1/2 


AD-A162  518  INCORPORATION  OF  ACTIVE  ELEMENTS  INTO  THE  ARTICULATED 
TOTAL  BODV  MODEL  <U>  PENNSVLVANIA  STATE  UNIV 
UNIVERSITV  PARK  DEPT  OF  INDUSTRIAL  AN  A  FREIVALDS 
UNCLASSIFIED  S8  JUN  85  AAMRL-TR-85-861  F22815-82-C-8586  F/G  fi/16 


NL 


FfT7-  fi 

„7j  jrJ  ;  .  *  ■'  <  ?  kfl 

hM 

?».r/.T >1  ■ '  *r7’  m 

r::  &  i  'Li”  m 

1  1  •  .  • 

k  $ 

JjT*  -rfRJK'i.  cWWf® 

gp1  f  j  1  4*  1  .  V'l.  1 

Vi  j 

>  7j i 

V.*  .;  -•  ..  _  •;-  ;  '  ■.  ..  p  .  .7).  . 

•-•’•  '•; i”# ''2fSi 

•;  .".  . .  s'  .  i,*i'  :  ,  . .  -v .  •  .. : :  ,  ■  •••'*5:  .  :  "  «  . 

§^®s£*”fcj’5 ' :' rf'rissh s'|"i*5f 


vv  -  ■?;V.t'*'^\-:^:  * ~-r 

f  '•  k  «i:-  •  '  '  t  .  i\  • 


HENNING  ESAON  GIERKE,  Dr  Ing 
Director 

Biodynamics  and  Bioengineering  Division 

Air  Force  Aerospace  Medical  Research  Laboratory 


gMI 

-i  ;  $ 

fTT  C^T^TTv’-  [  ••  - 

1  i?)  j  4 

iW? :  -_1 " 

WWfcA,  ’■*" 

{e®*  swT  v- 

vawt*£«»  « 

'  •••  •*•  ■  /• 

i*  V 

jcV^a 

f|i 

KPy?.^  0 

ja-y-'V  y*f  O^tr 

$*K$|g 

W\ 

j£Sj]| 

WfimSm 

SECURITY  classification  of  this  page 


REPORT  DOCUMENTATION  PAGE 


Is  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2».  security  classification  authority 


2b.  OECLASSI f ication/oowng  raoing  schedule 


3.  distribution/a vailability  of  report 
Approved  for  public  release. 
Distribution  Unlimited. 


4  PERFORMING  ORGANIZATION  REPORT  NUM8ERISI 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


AAMRL -TR-8  5 -061 _ 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

The  Pennsylvania  State 
University 


6c.  AOORESS  ICily.  Stale  and  ZIP  Code) 

Department  of  Industrial  &  Management 
Systems  Engineering 
University  Park  PA  16802 


b.  OFFICE  SYMBOL 
(If  applicable) 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Harry  G.  Armstrong  Aerospace  Medical 
Research  Laboratory 


7b.  ADDRESS  (City.  Stale  and  ZIP  Code) 

Wright-Patterson  AFB  OH  45433 


8a.  NAME  OF  FUNOING/SPONSORING 

8b.  OFFICE  SYMBOL 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION 

NUMBER 

ORGANIZATION 

(If  applicable) 

US  Air  Force/AFSC 

AAMRL 

F33615-83-C-0506 

8c.  ADDRESS  (City.  State  and  ZIP  Code) 

|  10.  SOURCE  OF  FUNOING  NOS. 

Aeronautical  Systems  Division 

PROGRAM 
ELEMENT  NO. 

PROJECT 

NO. 

TASK 

NO. 

WORK  UNIT 
NO. 

Wright-Patterson  AFB,  Ohio  45433 

61102F 

2312 

V3 

111.  TITLE  (Include  Security  Classification) 

62202F 

7231 

23 

07 

See  rpverse  side 

12.  PERSONAL  AUTHORIS) 

Freivalds,  Andris 


13*.  TYPE  OF  REPORT 

Final 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVEREO 

from  83Sep26  to  85Jun30 


14.  DATE  OF  REPORT  (Ye.  Mo..  Day)  I  15.  PAGE  COUNT 

1985  June  30  97 


17. 


FIELD 


06 


COSAT  I  CODES 


GROUP 


2  _ 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

Articulated  Total  Body  Model,  Biodynamics,  Biokinematic 
Modeling,  Neuromusculature 


19.  ABSTRACT  (Continue  on  reverae  if  neceaaary  and  identify  by  block  number > 

The  Articulated  Total  Body  (ATB)  Model,  based  on  rigid-body  dynamics  with  Euler  equations 
of  motion  and  Lagrange  type  constraints,  is  used  by  the  Harry  G.  Armstrong  Aerospace  Medical 
Research  Laboratory  to  predict  the  forces  and  motions  experienced  by  air  crew  personnel  in 
typical  flight  operations.  To  provide  a  more  realistic  representation  of  human  dynamics, 
active  neuromusculature  was  added  to  die  ATB  Model  using  elements  of  the  newly  developed 
advanced  harness  system, 

A  lumped  three-parameter  muscle  model  with  a  contractile  element,  a  damping  element 
and  a  parallel  elastic  element  was  developed.  The  contractile  element  included  a  length- 

tension  relationship,  a  force-velocity  relationship  and  an  active  state  function.  The 
basic  fiber  mechanisms  were  integrated  into  muscle  systems  utilizing  motor  unit  organization, 
orderly  recruitment  of  motor  units  and  adjustments  in  force  due  to  fatigue  and  reflex 
action.  The  complete  muscle  systems  were  then  used  to  replicate  the  human  neuromusculature 
of  the  trunk  and  neck  and  for  the  elbow,  shoulder,  hip  and  knee  joints*. 


20.  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT 

UNCLASSI  FIE  D/UNLIMITE  D  B  SAME  AS  RPT.  □  OTIC  USERS  □ 


22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

CORY  CARROLL 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


Unclassified 


22b  TELEPHONE  NUMBER 
(Include  Area  Code) 

(513)  255-5963 


22c  OFFICE  SYMBOL 

AAMRL/BBM 


SECURITY  CLASSIFICATION  OF  THIS  FAOE 


Block  11  continued. 

Incorporation  of  Active  Elements  into  the  Articulated  Total  Body  Model 


Block  19  continued. 


Several  validation  studies  were  performed.  One  simulated  elbow  flexion  with  one 
muscle/harness  system  representing  the  biceps  brachii  and  the  brachialis.  The  results 
indicated  that  the  force  velocity  effects  produced  the  greatest  changes  in  force,  with 
significant  force  changes  due  to  the  parallel-elastic  element.  Other  studies  simulated 
the  whole  body  response  to  a  2Gy  lateral  force  utilizing  trunk  musculature.  Although 
the  musculature  did  not  completely  prevent  the  lateral  deflection  of  the  body,  the 
response  was  significantly  delayed  compared  to  a  control  response,  with  the  head  and 
neck  maintaining  the  upright  posture  for  a  longer  period  of  time. 


SECURITY  CLASSIFICATION  OF  THIS  FAOE 


PREFACE 


Research  in  this  final  report  was  performed  under  United 
States  Air  Force  Research  Contract  No.  F33615-83-C-0506  awarded 
to  The  Pennsylvania  State  University  from  September  26,  1983 
to  June  30,  1985.  The  United  States  Air  Force  Project  Engineer 

was  Captain  Cory  D.  Carroll  of  the  Modeling  and  Analysis  Branch, 
Biodynamics  and  Bioengineering  Division,  Harry  G.  Armstrong 
Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB ,  Ohio 
The  author  expresses  his  sincere  thanks  to  Dr.  Ints  Kaleps  of 
the  Harry  G.  Armstrong  Aerospace  Medical  Research  Laboratory,  for 
his  suggestions,  evaluations  and  impetus  for  the  research. 


TABLE  OF  CONTENTS 


Page 

LIST  OF  FIGURES .  vi 

LIST  OF  TABLES . viii 

LIST  OF  SYMBOLS .  ix 

I.  INTRODUCTION  .  1 

II.  OBJECTIVES .  2 

III.  PHASE  I  -  THE  BASIC  MUSCLE  MODEL .  3 

A.  Skeletal  Muscle  .  3 

B.  Passive  Viscoelastic  Elements .  5 

C.  Active  Muscle  Elements  .  11 

IV.  PHASE  II  -  ADVANCED  DEVELOPMENT .  19 

A.  Organization  of  Fibers  into  Motor  Units  .  19 

B.  Orderly  Recruitment  Patterns  .  21 

C.  Feedback  Control  .  26 

D.  Time  Varying  Effects .  29 

V.  PHASE  III  -  MODELLING  THE  GENERAL  MUSCULATURE .  31 


A.  Elbow  Joint .  31 

B.  Shoulder  Joint .  33 

C.  Hip  and  Knee  Joint  .  . .  34 

D.  Trunk  and  Neck  Musculature .  34 

VI.  PHASE  IV  -  SOFTWARE  IMPLEMENTATION  OF  THE  NEUROMUSCULATURE  ....  35 

A.  Increasing  the  Number  of  Muscles .  35 

B.  Active  Neuromusculature  .  39 

VII.  SIMULATION  AND  VALIDATION .  44 

A.  Simulation  of  Trunk  and  Neck  Musculature .  44 

B.  Sensitivity  of  Muscle  Origins  and  Insertions  .  45 

C.  Anatomical  Studies  of  Trunk  Musculature  .  54 

VIII.  CONCLUSIONS .  54 

APPENDIX  A  -  INPUT  DATA  FOR  HUMAN  MUSCULATURE . 

Table  1  -  Elbow  Musculature . A-l 

Table  2  -  Shoulder  Musculature . A-2 

Table  3  ~  Hip  and  Knee  Musculature . A-3  -  A-4 

Table  4  -  Neck  Musculature  . A-5  -  A-6 

Table  5  -  Trunk  Musculature . A_7  _  a-8 


APPENDIX  B  -  LISTING  OF  FORTRAN  IV  SOURCE  DECK  FOR  SUBROUTINE 


MUSCLE .  B- 

APPENDIX  C  -  LISTING  OF  FORTRAN  IV  SOURCE  DECK  FOR  SUBROUTINE 

HINPUT .  C- 

APPENDIX  D  -  REVISIONS  TO  VALIDATION  OF  THE  CRASH  VICTIM  SIMULATOR 

VOL.  3,  USER'S  MANUAL . D- 

APPENDIX  E  -  ATB  SPECIFICATIONS  FOR  HUMAN  MUSCULATURE  . 


Table  6  -  Elbow  Musculature . ’if?- 

Table  7  -  Shoulder  Musculature . E- 

Table  8  -  Hip  and  Knee  Musculature . E_ 

Table  9  -  Neck  Musculature . E- 

Table  10  -  Trunk  Musculature . E- 

^‘APPENDllX  F  -  REFERENCES . . F- 


J^ccesion  7oT 

NT,S  CRa&T 

OTIC  TAB 

U-iannou  ,Ced 
Justification 


By . 

t  ib.  tio  ./ 


*wai,ato%Codes 

I  a  ri  I  - 


LIST  OF  FIGURES 


gure 


1  Schematic  Arrangement  of  Skeletal  Muscle  Fiber  Arrangement 


2  Molecular  Substructure  of  Mammalian  Skeletal  Muscle  .  .  . 


3  Lumped  Model  of  Skeletal  Muscle 


4  Simplified  Muscle  Model 


5  Passive  Muscle  Force-Strain  Function 


6  Passive  Muscle  Viscous-Damping  Forces 


7  Muscle  Length-Force  Relationship 


8  Muscle  Force-Velocity  Relationship 


9  Relative  Force  Developed  by  Active  State  Function  .  .  .  . 


10  Force  Developed  by  Orderly  Recruitment  of  Motor  Units  .  . 


1 1  Four  Types  of  Feedback  Control 

a)  Stretch  reflex  . 


b)  Reciprocal  inhibition  .  28 

c)  Clasp  Knife  Reflex  .  28 

d)  Flexor  and  crossed  Extensor  Reflex . 28 

12  Endurance  Time  as  a  Function  of  Relative  Strength . 32 


13  Flowchart  of  Subroutine  MUSCLE 


14  Graphical  Display  of  Trunk  Response  With  Advanced  Neuromusculature.  46 


15  Graphical  Display  of  Trunk  Response  to  Lateral  G-Forces  .  .  47 


16  Angular  Displacement  of  Upper  Trunk  Due  to  Lateral  G-^orces  ....  46 


17  Anatomy  of  the  Elbow  Flexors . 50 


18  Force  Production  by  the  Brachioradialis  during  Eccentive  Motion  .  51 


19  Force  Production  by  the  Brachial is  during  Eccentric  Motion  .  52 


20  Force  Production  by  the  Biceps  during  Eccentric  Motion 


vii 


LIST  OF  TABLES 


Specifications  on  Elbow  Musculature 


Specifications  on  Shoulder  Musculature 


Specifications  on  Hip  and  Knee  Musculature  .  .  . 


Specifications  on  Neck  Musculature 


Specifications  on  Trunk  Musculature 


ATB  Specifications  for  Shoulder  Musculature.  .  . 


ATB  Specifications  for  Shoulder  Musculature.  .  . 


ATB  Specifications  for  Hip  and  Knee  Musculature. 


ATB  Specifications  for  Neck  Musculature 


ATB  Specifications  for  Trunk  Musculature  .  . 


LIST  OF  SYMBOLS 


cross-sectional  area  of  muscle 

viscous  damping  coefficient 

cross-bridge  elastic  element 

contractile  element 

damping  element 

viscoelastic  spring  constant 

force  developed  in  muscle  normalized  by  the  maximum  Isometric 
tension 

force  developed  by  damping  element  normalized  by  the  maximum 
isometric  tension 

force  due  to  length-force  relationship  normalized  by  the  maximum 
isometric  tension 

force  developed  by  parallel  elastic  element  normalized  by  the 
maximum  Isometric  tension 


force  due  to  force-velocity  relationship  normalized  by  the 
maximum  isometric  tension 

force  developed  in  muscle 

maximum  isometric  tension 

spring  constant  for  parallel  elastic  element 

spring  constant  for  series  elastic  element 

instaneous  muscle  length 

resting  muscle  length 

mass 

relative  number  of  recruited  motor  units 

relative  number  of  Type  I  motor  units 

relative  number  of  Type  II  motor  units 

number  of  recruited  motor  units 

number  of  Type  I  motor  units 


number  of  Type  II  motor  units 


total  number  of  motor  units 


functional  form  for  muscle  model 


rate  of  change  of  active  state 


parallel  elastic  element 


parallel  elastic  element  of  the  3arcolemma 


active  state  function 


slack  in  muscle  length 


series  elastic  element 


time 


contracation  time  of  motor  unit 


endurance  time 


*f-T0 


resting  Ca++  concentration 


free  Ca++  concentration 


I  Yf 

e  »  muscle  strain 

I 

e  -  muscle  strain  rate 

^MAX  ■  maximum  muscle  strain  rate 

e'  -  pretensioned  muscle  strain 

n  -  viscous  damping  coefficient 

n  -  normalized  muscle  velocity 

X  -  length  of  contractile  element 

X  -  resting  length  of  contractile  element 

j  -  maximum  shortening  velocity  of  contractile  element 

*0 

\q  -  length  of  cross  bridge  elastic  element 

Xg  -  length  of  series  elastic  element 

u  -  microns 


relative  stimulation  rate 


normalized  length  of  contractile  element 


stimulation  rate 


maximum  stimulation  rate 


xiii 


I.  INTRODUCTION 


Biodynamic  computer-based  models  for  the  prediction  of  human  body 
response  to  mechanical  stress  have  become  extremely  useful  and  cost-effective 
research  and  developmental  tools,  especially  as  alternatives  to  direct 
experimentation  with  humans  and  animals.  These  models  attempt  to  simulate  or 
predict  the  forces  and  motions  experienced  by  a  body  in  high-acceleration 
events  such  as  impacts  or  from  sudden  forces  such  as  wind  shear.  In 
particular,  the  Air  Force  is  interested  in  the  reactions  of  aircrew  personnel 
to  such  forces  typically  encountered  in  various  phases  of  flight  operations, 
including  emergency  ejections  from  high-speed  aircraft.  Such  a  hazardous 
environment  is  well  suited  to  computer  modeling,  and  with  proper  execution, 
considerable  insight  into  body  motion  and  stresses  developed  in  the  body  can 
be  gained. 

The  Modelling  and  Analysis  Branch  of  the  Biodynamics  &  Bioengineering 
Division  of  the  Air  Force  Aerospace  Medical  Research  Laboratory  ( AFAMRL)  has 
been  using  a  human  body  modelling  computer  program  known  as  the  Articulated 
Total  Body  (ATB)  Model  for  several  years.  The  model  is  based  on  rigid-body 
dynamics  using  Euler  equations  of  motion  with  Lagrange-type  constraints  (Fleck 
et.al.  197*0.  The  specific  configuration  uses  15  body  segments  (head,  neck, 
upper  torso,  center  torso,  pelvis,  upper  arms,  lower  arms,  upper  legs,  j.ower 
legs,  and  feet)  and  14  joints  between  the  segments  (Fleck  and  Butler,  1975). 
Although  it  was  originally  developed  by  the  Calspan  Corporation  for  the  study 
of  human-body  and  anthropometric-dummy  dynamics  during  automobile  crashes  for 
the  United  States  Department  of  Transportation  (Fleck  et.al.  1974;  Fleck, 
1975),  the  ATB  Model  was  sufficiently  general  to  allow  simulation  of 
whole-body  articulated  motion  resulting  from  various  impacts  or  abrupt 
accelerations  applied  to  the  body.  Furthermore,  modifications  involving 
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special  joint  forces,  aerodynamic  forces  and  a  complex  harness  system  were 
added  to  accommodate  specific  Air  Force  applications  (Fleck  and  Butler,  1975). 

The  ATB  Model  initially  reflected  human  body  structure,  mass  distribution 
and  tissue  material  properties  for  passive  responses.  An  early  effort  to 
improve  the  ATB  Model  in  regards  to  active  responses  resulted  in  the 
development  of  a  lumped  three  parameter  viscoelastic  muscle  model  superimposed 
on  the  advance  restraint  system.  (Freivalds,  1981!;  Freivalds  and  Kaleps, 

1983;  Freivalds  and  Kaleps,  198^) .  However,  the  early  efforts  were 
constrained  by  the  low  number  (five)  of  harness  systems  provided  in  the  ATB 
Model,  limiting  simulations  to  simple  joint  motions  or  very  crude  whole  body 
motion.  Also,  complex  neuromuscular  functions  such  as  motor  unit  recruitment 
patterns,  time  varying  effects,  etc.  were  not  included.  Thus,  further 
development  of  the  neuromuscular  system  was  needed  to  better  simulate  active 
human  responses  to  high  G  forces. 

II.  OBJECTIVES 

The  objective  of  this  project  was  to  further  define  and  formulate 
methodologies  for  implementing  active  muscle  responses  into  the  present  ATB 
Model.  Three  considerations  were  Involved:  (1)  basic  muscle  phenomena  such 
as  motor  units,  recruitment  patterns,  and  fatlque  were  to  be  included,  (2) 
particular  emphasis  was  to  be  placed  on  muscles  acting  in  the  torso  and  neck 
region  which  affect  flexion,  extension  and  lateral  motion  of  the  trunk  in  a 
seated  posture  and  (3)  the  modifications  should  be  transparent  to  the  user, 
such  that  existing  input  decks  would  remain  valid. 

The  objective  was  approached  in  four-phased  approach.  In  Phase  I,  the 
basic  muscle  model  developed  during  the  early  efforts  (Freivalds,  1982; 
Freivalds,  1 984 ) ,  was  re-examined  and  redefined.  In  Phase  II,  advanced 
features  of  types  of  motor  unit,  motor  unit  recruitment,  force  buildup  and 
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endurance  times  were  developed  and  included  into  the  ATB  Model.  In  Phase  III, 
the  ATB  Model  was  modified  to  allow  for  up  to  50  muscles  and  a  representative 
musculature  for  the  entire  body  was  developed.  In  Phase  IV,  various 
simulations  were  performed  in  order  to  validate  the  modelling  efforts. 

III.  BACKGROUND 

A.  Skeletal  Muscle: 

Skeletal  muscles  usually  originate  on  the  skeleton,  3pan  one  or  more 
Joints  and  insert  into  a  part  of  the  skeleton  again.  Each  muscle  is  enclosed 
in  a  connective  tissue  sheath  called  the  epimysium  and  is  held  in  its  correct 
position  in  the  body  by  layers  of  fascia.  The  muscle  is  attached  to  the  bones 
via  tendons,  while  the  interior  is  compartmentalized  into  longitudinal 
sections  called  the  fasciculi,  each  containing  many  individual  muscle  fibers. 
The  fibers  are  enveloped  by  a  connective  tissue  called  the  endomysium,  which 
transmits  the  force  of  the  muscle  contraction  from  individual  fibers  to  the 
tendons  (Fung,  1981). 

The  muscle  fibers  do  not  always  run  parallel  to  the  force  transmitting 
tendons,  as  they  do  in  fusiform  muscles.  They  can  be  arranged  in  unipennate, 
bipennate  or  multipennate  form,  thus  altering  the  force  transmittinig 
characteristics  (Fig.  1). 

The  muscle  fiber,  the  basic  structural  unit,  with  a  diameter  of  10-60  p 
and  length  from  several  millimeters  to  several  centimeters,  can  be  subdivided 
further  into  myofibrils  of  Ip  diameter.  These  myofibrils  comprise  the 
hexagonal  array  of  protein  filaments  that  are  directly  reponsible  for 
the  contractile  process  and  give  rise,  with  appropriate  stains  to  the  peculiar 
striations  that  are  characteristic  of  skeletal  muscle  (Figure  2).  A  repeating 
unit  known  as  the  sarcomere  is  defined  by  the  vertical  z-disk.  Two  types  of 
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8IPENNATE 


MULTI  PENNATE 


Fig.  1 


Fig.  2 


Schematic  representation  of  skeletal 
muscle  fibre  arrangement. 


z-disk 


sarcomere 


Molecular  substructure  of  mammalian 
skeletal  muscle. 


4 


protein  filaments  are  distinguishable  in  each  sarcomere,  thin  ones  about  5nm 
(50A)  in  diameter  and  thicker  ones  about  12nm  (120A)  across.  The  thin 
filaments  contain  actin,  globular  molecules  in  a  triple  helix,  while  the  thick 
filaments  contain  myosin,  long  molecules  with  globular  heads.  The  thin 
filaments  are  each  attached  at  one  end  to  a  z-disk  and  are  free  at  the  other 
to  interlace  with  the  thick  filaments.  The  A-band  is  the  region  of  overlap 
between  thick  and  thin  filaments,  the  I-band  contains  solely  the  thin 
filaments,  while  the  H-band  is  the  middle  region  of  the  A-band  into  which  the 
actin  filaments  have  not  penetrated  (Fung,  1981). 

The  actual  contractile  process  takes  place  at  the  Junctions  between  the 
myosin  and  actin  in  a  process  known  as  the  sliding  filament  theory  first 
presented  by  H.  E.  Huxley  (1953).  The  myosin  molecules  consist  of  a  long  tail 
piece  and  a  "head".  The  tails  lie  parallel  in  a  bundle  to  form  the  core  of 
the  thick  filament  while  the  heads  project  laterally  from  the  filament  in 
pairs,  rotated  with  respect  to  its  neighbors  to  form  a  spiral  pattern  along 
the  filament.  These  heads  seem  to  be  able  to  nod;  they  lie  close  to  their 
parent  filament  in  relaxation,  but  stick  out  to  actin  filaments  when  excited. 
Thus,  during  muscle  contraction  the  muscle  fiber  shortens  as  the  filaments 
slide  over  each  other,  forming,  breaking  and  reforming  chemical  bonds  between 
the  myosin  heads  and  the  globular  actin  molecules. 

B.  Passive  Viscoelastic  elements: 

Previous  modelling  efforts  (Freivalds,  1982;  Freivalds,  1984;  Freivalds 
and  Kaleps,  1983;  Freivalds  and  Kaleps,  1984)  produced  the  lumped  model  of 
skeletal  muscle  shown  in  Figure  3.  Structures  which  lie  in  parallel  to  the 
force  producing  sarcomeres:  the  sarcolemna  (sheath)  of  the  individual  fiber 
and  the  various  outer  connective  sheaths  (fascia,  endomysia,  perimysia)  are 
represented  by  the  parallel  elastic  element  (PE).  Practically  all  the  tension 


observed  when  stretching  the  resting  muscle  will  result  from  this  element. 
Because  the  muscle  fiber  is  60$  water,  an  appropriate  damping  element  (DE)  is 
also  included  parallel  to  the  contractile  elements  (CE).  The  contractile 
element  represents  the  purely  contractile  protein  molecules.  In  series  with 
the  contractile  element  is  a  bridge  element  (BE)  representing  the  elastic 
elements  within  the  cross  bridges  and  the  z-disks.  The  parallel  elastic 
element  for  the  sarcomere  (PS)  does  not  contain  a  damping  component,  since  the 
sarcolemma  attached  to  the  z  disks  does  not  allow  appreciable  movement.  The 
tendinous  parts  of  the  muscle  fiber  are  located  near  the  origin  and  insertion 
of  the  fiber  and  thus  are  depicted  by  a  series  elastic  element  (SE).  Any  mass 
of  the  sarcomeres  is  disregarded,  especially  when  compared  to  the  much  larger 
external  mass  that  the  muscle  contraction  must  move  (Hatze,  1981). 

Certain  assumptions  allow  a  further  reduction  of  the  lumped  model  to  a 
simpler  model  yet.  SE  and  BE  can  be  considered  to  be  very  stiff  springs  and 
eliminated  completely.  This  contention  is  supported  by  Bawa  et.al.  (1976)  who 
found  Ksg-3724  N/m  to  be  much  larger  than  Kpg-1000  N/m.  KBg  can  be  considered 
to  be  in  a  similar  range  with  Kgs*  Eliminating  SE  and  BE  results  in  a  model 
with  four  parallel  elements,  two  of  which  are  elastic  and  can  be  combined  into 
one  parallel  elastic  element.  The  final  simplified  model  is  given  in  Figure 
4. 

The  total  force  developed  by  the  simplified  model  of  Figure  4  can  now  be 
expressed  as: 

F  -  ( fpE  +  fcE  +  fDE)FMAX  (D 

where  F^x  I3  the  maximum  isometric  tension  of  the  muscle. 


For  concentric  or  shortening  contractions  both  fpg  and  fpg  zero,  PE 
producing  force  only  under  stretch  and  DE  producing  force  only  under  external 
tension,  i.e.  eccentric  motion. 

Mathematical  representations  for  each  element  were  determined  as  follows. 
For  the  parallel  elastic  element,  extensive  tests  on  the  tensile  properties  of 
resting  human  sartorius  muscle  carried  out  by  Yamada  (1970)  indicate  an 
exponential  force-strain  function: 

fPE  -  .ooi6296(e7-66l6e-i )  (2) 

where  fpE  is  the  force  developed  by  the  PE  normalized  with  respect  to  maximum 
isometric  tension  in  the  muscle  and  e  is  the  strain: 


where  l  is  the  instananeous  muscle  length  and  H0  is  the  resting  length.  This 
force-strain  curve  is  shown  in  Figure  5. 

The  velocity  dependence  of  the  damping  element  (DE)  can  be  expressed 
similarly  to  the  form  used  for  a  simple  mechanical  dashpot: 


fD£  -  n  e 


where  foe  the  normalized  force,  e  is  the  muscle  strain  rate  and 
n  is  the  viscous  damping  coefficient.  As  mentioned  previously  (Frelvalds, 
1984)  the  force  produced  by  the  passive  viscous  damping  element,  amounted  to 
as  much  as  a  30%  increase  above  MVC.  This  seemed  to  be  an  excessively  large 
effect  for  a  passive  response  and  consequently  the  damping  coefficient  was 
re-evaluated.  Alexander  and  Johnson  (1965)  loaded  the  frog  sartorius  muscle 
dynamically  in  both  the  passive  and  active  states.  From  the  ramp  shaped 
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loading  curves  one  can  establish  time  constants  for  the  resulting  strain 
responses  as  follows.  Using  viscoelastic  modelling  techniques  (Bland,  I960; 
Flugge,  1967)  one  can  express  the  relationship  of  a  parallel  combination  of  a 
spring  and  dashpot  as: 

o(t)«  Ee(t)  +  ne(t)  (5) 

where  o  -  stress,  and  E  ■  spring  constant.  Taking  the  Laplace  transform  of 
both  sides  yields: 

o(s)  -  (E  +  ns)  e(s)  (6) 

If  the  loading  function  is  a  ramp  function  of  the  form  OQt,  then  its  Laplace 
transform  is  oq/s2  and  the  resulting  strain  becomes: 


(7) 


(8) 


(9) 


Based  on  the  work  of  Alexander  and  Johnson  (1965),  the  strain  n  in  Eq.  9 
can  be  considered  a  composite  strain  consisting  of  a  dynamic  strain  and  a 
passive  strain.  At  time  t  equal  H  sec,  the  resulting  elongation  from  a 


linear  ramp  loading  of  4.4  gr/sec  is  .5  mm  for  the  passive  state  and  .075  mm 
for  the  active  state.  In  the  passive  state  the  viscous  damping  forces  for 
both  the  contractile  element  and  secondary  tissues  come  into  play.  During  the 
active  state  the  damping  effect  in  the  contractile  element  is  overshadowed  by 
the  force  velocity  relationship  and  only  the  damping  effect  due  to  other 
tissues  is  measured.  Based  on  the  data  of  the  frog  sartorius  E  «  6.57  gr/sec, 
oq  -  4.4  gr/sec,  t  -  4  sec).  Alexander  and  Johnson  (1965)  found  a  value  of 
6.38  for  the  ratio  of  the  combined  damping  coefficient  (found  during 
the  passive  state)  to  the  passive  damping  coefficient  (found  during  active 
state) . 

Converting  the  cat  plantaris  combined  damping  coefficient  of  63  Nsec/m 
(which  is  more  similar  to  human  values  as  opposed  to  the  frog  coefficient) 
found  by  Bawa,  Manna rd  and  Stein  (1976)  yields  a  value  of  9.87  Nsec/m. 
Normalizing  to  the  plantaris  length  of  .05  m  and  maximum  contraction  force  of 
2.45  N  yields  the  final  passive  damping  coefficient  value  of  .139  sec.  The 
complete  functional  relationship  of  Eq.  4  is  plotted  in  Fig.  6.  Note  that 
damping  is  only  effective  on  eccentric  movements. 

C.  Active  Muscle  Elements: 

The  contractile  element  is  the  only  active  component  in  the  model.  Its 
behavior  is  extremely  complex  and  depends  nonlinearly  on  its  length, 
contractive  history,  velocity  of  movement,  the  degree  of  stimulation  and  its 
temperature.  However,  for  practical  purposes,  only  three  basic  functions  were 
considered:  the  length-force  relationship  (f^),  the  force-velocity 
relationship  (fv),  and  the  active  state  function  (fq).  Thus 
fee  *  fq  •  •  fv 
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The  length-force  relationship  is  determined  by  the  number  of  active  cross 
links  or  filamentary  overlap  and  can  be  adequately  expressed  from  the  data  of 
Gordon  et.al.  (1966)  by  the  function  suggested  by  Hatze  (1981,  p.  42): 

fj,(e)  -  .32  +  .71e_1-112e  sln(3.722(e+.344))  (1  1  ) 

This  function  is  shown  in  Figure  7  and  further  details  are  given  in  Freivalds 
(1982). 

The  force-velocity  relationship  is  determined  by  the  rate  of  breaking  and 
reforming  the  cross  bridges  with  higher  rates  producing  less  effective  bonds. 
To  account  for  the  whole  range  of  negative  velocities  (shortening  or 
concentric  contractions)  as  well  as  positive  velocities  (lengthening  or 
eccentric  contractions)  Hatze  (1981,  p.  45-46)  has  defined  the  following 
expression: 

f  ( n)  -  .1433  1*1073  ♦  e"1'1409  3inh  (3.2n  +  1.6)  }_1  -.005[2-e6e]  (12) 


where  fv(n)  is  the  normalized  force  due  to  the  force-velocity  relationship 
as  defined  by  the  first  term  and  reduced  by  internal  resistance  as  defined  by 
the  second  term.  However,  since  the  coefficient  of  the  second  term  is  smaller 
by  a  factor  of  30  than  the  first  term,  it  can  be  disregarded  for  present 
purposes,  n  represents  the  normalized  contractile  element 
velocity: 


with  being  the  maximum  shortening  velocity  of  the  contractile  element. 

Equation  6  represented  by  Figure  8,  with  further  details  given  in  Freivalds 


CONCENTRIC  ECCENTRIC 


The  active  state  function  fq  is  defined  by  the  relative  amount  of  calcium 
(Ca*+)  bound  to  troponin  (inhibitor  molecule  of  actin).  If  the  maximum 
number  of  potential  interactive  sites  on  the  thin  filament  are  exposed  by  the 
action  of  Ca++,  then  q-1 ;  while  in  a  resting  state  q-q0.  Thus  the  Isometric 
tension  developed  by  a  muscle  fiber  at  a  given  length  £q  the  CE  is  directly 
proportional  to  q  (Hatze,  1981,  p.  33). 

Define  Y  to  be  the  difference  between  the  real  free  Ca++  concentration  Yf 
and  the  free  Ca++  concentration  Y0  in  the  resting  fiber.  However,  for 
practical  purposes  since  Y0  <<  Yf,  we  have  Y-Yf.  Let  p  -  dq/dY  denote  the 
Ca*+  concentration  rate  of  change  of  the  active  state  q.  The  process  of 
binding  Ca+*  ions  to  the  troponin  sites  is  hypothesized  by  Hatze  0981)  and 
supported  by  the  experimental  studies  of  Ebashi  and  Endo  (1968)  to  be  a 
function  of  the  length  8,  of  the  CE  and  of  the  difference  between  the  maximum 
and  present  value  of  q  and  controlled  by  a  negative  feedback  loop  as 
follows. 

dp/dY  -  pi2(e)  (1-q)-2p2Pi (e)p  (14) 

where  e  -  ( IL-JL0 ) / 8.0  is  the  strain  and 

dq/dY  *  p  (15) 

Solving  the  differential  system  of  Equation  14  and  Equation  15  with  initial 
conditions: 
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* - N 
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(16) 
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(17) 

One  obtains  a  normalized  solution: 


(1-q0> 

V  '  "V 


.  ra2p1(e)Y  m^CejY 

i 


where 


/  2  ,J/2  >.  , 

mi ,2  ”  ~p2  -  p2  ~  P2  >  1 


(II 


(1  ‘ 


Substituting  experimentally  found  values  (Hatze,  1981): 


2.3**  x  10 
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(£ULi) 

v  e+1  ’ 


1/2 


-.44  <  e  <  .8 


(2( 


p2  *  1.05 


(2' 


One  obtains: 


q(e,Y) 

I 


-  (,-,0)(2.1»  e*’-167  -  ’°T  hU)m) 

_1  l4e-2.096  x  107h(e)Y(t) ) 


(2; 


where 


(2: 


h(  e) 


( 


£■*■.44. 

e+.r 


1/2 


However,  a  simple  computational  approximation  is  provided  by  Hatze  (1981,  p 
40)  for  most  mammalian  muscles: 


2  2 

q  +P 

q(e,Y)  -  - r - —  (2 

1+P  (ejY^ 

I 


where 


p(e)  -  66,2' 


1 . 9(  £•*•1 ) 


The  function  Y,  the  free  Ca+  +  ion  concentration  can  be  represented  as  a 
function  of  time  t  and  the  stimulation  rate  v  by  a  trend  function  which 
represents  the  average  behavior  of  Y  in  successive  time  intervals,  and  which 
approaches  a  maximum  value  asymptotically  and  has  the  rate  of  increase 
proportional  to  the  stimulation  rate  (Hatze,  1981,  p.39). 

Y  »  m(Cv-Y)  Y(0)  -  Y0  (26) 

where  m  and  C  are  constants  (C  -  1.373  x  1C1*  (Hatze,  1981)  and  m  to  be 
determined  later)  and  v  is  the  relative  stimulation  rate  defined  by 

0  <  v  -  ^  <1  (27) 

where  t”1  and  t-1  denote  the  stimulation  rate  and  maximum 
stimulation  rate  respectively. 

The  rate  of  stimulation  of  motor  units  during  voluntary  contraction  has 
been  very  controversial.  Several  studies  have  found  a  fairly  constant 
discharge  frequency  over  a  wide  range  of  tension  for  individual  motor  units 
(Bigland  and  Lippold,  1 954 ;  Clamann,  1970),  while  others  maintain  that  an 
Increase  in  muscle  tension  is  achieved  in  part  by  an  increase  in  the 
stimulation  rate  and  that  this  may  be  important  in  achieving  precision  and 
smoothness  of  contraction  (Marsden  et.al  1971;  Person  and  Kudina,  1972; 
Milner-Brown  et.al  1973b).  However,  even  for  those  studies  who  found  a  rise 
in  discharge  frequency  with  tension,  the  frequency  at  the  start  of  the 
discharge  for  rapid  contractions  was  much  higher  and  closer  to  the  maximum 
stimulation  rate  (Tanji  and  Kato,  1973).  Thus,  it  is  fairly  reasonable  to 
assume  a  constant  stimulation  rate  and,  therefore,  a  constant  relative 
stimulation  rate. 


Now  solving  Equation  26  with  v  constant  yields: 


Y  -  (YQ-Cv)  e"mt  +  Cv 


With  Y0  <<  Cv  (Y0  -  1x10-9;  Hatze,  1981,  p.  33),  Equation  28  reduces  to: 


Y  -  Cv(1-e-mt) 


Substituting  Equation  29  and  Equation  25  (with  e  -  0)  into  Equation  24  yields  the 
active  state  function: 


q(t)  - 


.005  +  82.63  v2  (1 -e~mt) 
1  ♦  82.63  v2  (1-e~mt) 


which  is  represented  in  Fig.  9.  Consequently,  f^g  can  be  redefined  by  using 
the  relative  force  fq  developed  by  the  active  state  function  q: 


fCE  “  fq(t)  f*(e)  fv(ri) 


IV.  PHASE  II  -  ADVANCED  DEVELOPMENT 


A.  Organization  of  Fibers  into  Motor  Units 

The  population  of  motor  units  can  be  subdivided  into  two  distinct 
populations  based  on  their  contractile  and  histochemical  properties:  Type  I 
(slow  twitch)  motor  units  and  Type  II  (fast  twitch)  motor  units  (Close,  1972). 
Type  I  motor  units  have  slower  contraction  times,  tend  to  be  more  aerobic  and 
less  fatigable  and  are  recruited  at  lower  tension  levels.  Type  II  motor  units 
have  faster  contraction  times,  tend  to  be  more  anaerobic  and  more  fatigable 
and  are  recruited  at  higher  tension  levels  (Close,  1972;  Milner  -  Brown, 
et.al.,  1973a).  Thus,  the  total  muscle  force  output  should  be  the  sum  of  the 
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TYPE  n  M.U..REL.  STIM 


force  output  from  Nj  of  the  Type  I  motor  units  and  Njj  of  Type  II  motor 
units: 

F«(f  +  f  +  f  +f)F 

PE  CEj  CEjj  DE'  MAX  (32) 

Furthermore,  the  same  total  population  can  be  subdivided  into  two-dynamically 
different  populations:  the  N  population  of  active  motor  units  and  the  N-N 
population  of  inactive  or  resting  motor  units,  where  N  is  the  total  number 
of  motor  units  in  the  muscle. 

Muscle  properties  dependent  on  fiber  type  will  be  developed  in  the  next 
section  as  many  of  these  also  depend  on  the  recruitment  pattern  used. 

B.  Orderly  Recruitment  of  Motor  Units 

It  has  been  well  established  that  motor  units  are  recruited  in  a 
sequential  order  according  to  their  sizes  (Milner  -  Brown,  et.al.  1973a).  The 
cumulative  relative  cross-sectional  area  u  occupied  by  the  fibers  of  the 
recruited  units  increases  by: 

c  (N-1) 

u  «  u  e  — — -  o  <  u  <  u  <  1  (33) 

o  -  o 

where  uq  is  the  cross-sectional  area  of  the  smallest  motor  unit, 

c  =  -in  uQ,  (34) 

N  is  the  number  of  stimulated  motor  units  and  N  is  the  total  number 
of  motor  units  (Hatze,  1979).  For  N  large,  Equation  24  reduces 


u  e 
o 


(35) 


cN/N 


(1 


-N/N) 


(1-n) 


where  n  is  the  normalized  number  of  recruited  or  active  motor  units  «  N/N. 

The  Aj^u  of  the  relative  cross-sectional  area  u  upon  recruitment  of  the 
ith  motor  unit  is  then  defined  by: 


Atu  -  C  eCi/N  (36) 

where  C  is  a  normalization  constant  determined  by  the  requirement  that 
N 

l  Aiu  »  1  i.e. 
i-1 

C  .  ■■■_■ l-._  ■■  (37) 

N  c.  /N 
I  e1 
1-1 

Applying  the  ratio  of  the  smallest  to  the  largest  motor  cross-sectional  areas 
measured  in  a  muscle,  an  estimate  of  the  value  of  u0  for  a  given  muscle  can  be 
found.  These  range  from  uo».005  for  the  human  rectus  femorus  muscle  to 
uo*.009  for  the  human  biceps  muscle  (Hatze,  1979)  with  an  average  value  of 
u0=. 00673  resulting  in  c-5  to  be  used  for  the  present  study. 

Thus  two  very  important  properties  of  motor  unit  recruitment  dynamics 
have  been  included;  motor  units  are  normally  recruited  sequentially  from  the 
smallest  to  the  largest  and  the  size  of  the  recruited  units  as  well  as  the 
total  force  produced  grows  exponentially  (Fig.  10).  Combining  the  two  sets  of 
overlapping  population  distribution  of  motor  units  yields  two  distinct  cases: 
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c  - 
K- 


a)  N  <  Nj,  ie.  only  part  of  the  Type  I  motor  units  are  stimulated  and  none  of 
the  Type  II  can  be  stimulated  because  of  the  orderly  recruitment  pattern  and 

b)  N  >  Nj,  ie.,  all  Type  I  motor  units  are  stimulated  and  some  of  the  Type  II 
motor  units  are  stimulated,  but  none  of  the  Type  I  are  inactive.  These 
conditions  are  included  for  the  muscle  properties  for  which  they  are 
appropriate. 

Appropriate  properties  of  motor  units  can  be  obtained  from  the  data 
presented  in  Henneman  and  Olson  (1965).  The  contraction  time  tc  of  a  motor 
unit  is  a  decreasing  function  of  the  fraction  n  of  recruited  motor  units: 

tc  -  a2~a3n  (38) 

Thus  for  Type  I  motor  units 

\  “  a2I“a3In  0  <  n  <  nr  (39) 

and  for  Type  II  motor  units 

t  =  a„  -  a_  n  nT  <  n  <  1  (40) 

CII  2II  3II  I  -  - 

The  constants  a_  ,  a,  ,  a  ,  a.  can  be  determined  from  experimental 
dl  3I  2II  JII 


values.  For  n=*0,  the  value  of  tc  corresponds  to  the  contraction  time  of  the 
slowest  Type  I  unit  in  the  muscle,  approximatley  equal  to  .1  sec.;  for  n-nj 
the  value  of  tc  corresponds  to  the  fastest  Type  I  unit,  approximately  .045 
sec.;  and  for  n-1  (given  ni*1)  the  value  of  tc  corresponds  to  the  fastest  Type 
II  unit,  approximately  .01  sec.  (Grimby  and  Hannerz,  1977).  Substituting  and 
solving  for  the  unknown  values  yields: 
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.1 


I 


(41 ) 


a,  -  .055/nT 


(42) 


.01 


+ 


■  035 
(1-n:) 


a 


3 


II 


.035 

(1-nj) 


(43) 

(44) 


Substituting  Equations  41  and  42  into  Equation  39  and  Equations  43  and  44  into 
Equation  40  yield  respectively: 


,1  - 


.055  n 


0  <  n  <  n 


(45) 


'II 


.01  ♦  -035  iia 

1-n 


I 


nj  <  n  <  1 


(46) 


Several  other  important  parameters  can  be  derived  using  Equations  45  and 
46. Close  (1965)  showed  that  for  mammalian  skeletal  muscle  the  maximum  normalized 
speed  of  shortening  is  related  to  the  contraction  time  of  a  muscle,  consisting 
predominantly  of  one  fiber  type,  by: 


where  B  has  a  value  of  .297  for  human  muscle.  Using  Equations  45,  46  and  47, 


is  found  to  be  2.97/sec.  for  slow  (n-0)  and  29.7/sec.  for  fast 
(n-1)  motor  units.  However,  for  present  modelling  purposes,  an  average 
value  for  t-MAX  be  used.  Integrating  tc  from  Equation  45 

over  the  pattern  of  motor  unit  recruitment  for  Type  I  motor  units  (with  c  - 
5  for  human  muscle): 

(.1  -  .055x)  cx 

-  e  dx 

nI 

- i -  (48) 
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yields: 


e5n(5n-1 )-1 


(49) 


A  similar  process  of  integrating  over  the  pattern  of  Type  II  motor 
unit  recruitment  yields  an  average  contraction  time  of 


rr°n\  (e5n  (5n-1)  -  e  5(5n  -1))  (50) 

c  -  oi  ♦  ^235  _  (J__V _ : _ 

11  '~*i  (e5n-e5nl) 

The  rate  constant  In  Equation  29  also  depends  on  the  contraction 

time: 

m  -  A0/tc  (51) 

where  A0  -  .372  for  human  muscle  (Hatze,  1981,  p.  62).  For  slow  motor  units 
(n-0)  m  -  3.72  and  for  fast  motor  units  (n-1)  m  -  37.2.  Again  an  average 
value  is  used  for  present  modelling  efforts. 

C.  Feedback  Control 

Another  property  of  mammalian  skeletal  muscles  involves  the  feedback 
control  between  closely  related  muscles.  For  example  the  flexor  and  extensor 
muscles,  namely  the  biceps  and  triceps,  determine  in  unison  the  degree  of 
elbow  flexion.  The  afferent  neurous,  termed  la,  pass  information  from  the 
sensory  receptors  (usually  muscle  spindles)  in  the  muscle  to  the  spinal  cord 
where  they  form  monosynaptic  excitatory  connections  with  large  a  motoneurons 
which  carry  information  back  from  the  spinal  cord  to  muscle  fibers.  Thus 
stretch  of  one  muscle  will  cause  excitatory  impulses  in  the  spindle  to  excite 
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the  la  fiber  which  in  turn  exites  the  a-efferent  motoneurons  resulting  in  a 
reflex  contraction.  This  monosynaptic  reflex  arc,  shown  in  Fig.  11a,  is  known 
as  stretch  reflex.  However,  in  most  cases,  a  voluntary  active  response  will 
completely  override  in  magnitude  the  stretch  reflex.  Therefore  the  stretch 
reflex  hasn't  been  specifically  included  in  the  ATB  Model. 

A  second  reflex,  termed  reciprocal  inhibition,  involves  the  collaterals 
of  the  la  fibers  forming  disynaptic  inhibitory  connections  with  the 
antagonistic  motoneurons.  Activity  in  the  la  fibers  results  in  excitation  of 
one  set  of  motoneurons  (agonists)  and  at  the  same  time  inhibition  of  the 
antagonistic  motoneurons  as  shown  in  Fig.  11b.  Thus  the  reciprocal 
antagonistic  inhibition  facilitates  contraction  of  the  agonistic  muscle 
triggered  by  the  la  fiber  activity  by  simultaneously  inhibiting  the 
antogonistic  muscle  acting  on  the  same  joint  (Houk  and  Henneman,  1967; 

Schmidt,  1975;  Granit  and  Pompelano,  1979).  Typically  this  inhibition  appears 
to  be  fairly  complete  (Granit,  1970)  and  thus  a  reasonable  assumption  is  to 
eliminate  the  antagonist  muscles  completely  in  voluntary  motion  involving 
reciprocal  inhibition. 

It  is  also  possible  to  influence  the  output  of  the  muscle  spindle 
afferents  and  thus  adjust  the  amount  of  a-efferent  activity  returning  to  cause 
muscle  contraction,  by  a  second  set  of  efferent  nerves.  These  are  smaller  in 
size  termed  Y,  and  innervate  the  intrafusal  fibers  of  the  muscle  spindle. 
Increased  Y  activity  causes  an  increased  timing  frequency  of  the  la  fibers 
resulting  in  Increased  agonist  contraction.  Such  changes  represent  a  change 
in  the  gain  of  the  system  and  are  included  inherently  in  the  motor  unit 
recruitment  algorithm. 

Apart  from  the  muscle-spindle  receptors,  the  skeletal  muscle  contains 
another  type  of  receptor,  the  Golgi  tendon  organ,  that  is  important  for  a 


r 
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c.  CLASP  -  KNIFE  REFLEX 


d.  FLEXOR  ANO  CROSSEO 
EXTENSOR  REFLEX 


Fig.  11  Four  Types  of  Feedback  Control 


third  type  of  reflex.  The  Golgi  tendon  organs  are  located  at  the  tendons 
the  extrafusal  muscle  fibers  with  afferents  termed  lb  fibers  leading  to  the 
spinal  cord.  During  a  muscle  contraction  the  firing  frequency  of  lb  fibers 
increases  indicating  the  Golgi  tendon  organs  respond  to  tension.  In 
functional  terms  the  segmental  connections  of  the  a  lb  fibers  mirror  those  of 
the  la  fibers.  The  tendon  organs  form  inhibitory  connections  with  agonistic 
motoneurons  and  excitatory  connections  with  antagonistic  motoneurons  as  shown 
in  Fig.  11c.  A  strong  increase  in  muscle  tension,  whether  caused  by  passive 
stretch  or  active  contraction,  will  lead  to  inhibition  of  the  agonist 
motoneurons  via  lb  fibers  and  thus  prevent  a  too  strong  rise  in  tension. 

Under  severe  muscle  stretch,  muscle  tone  will  decline  suddenly,  resulting  in 
the  clasp-knife  reflex  which  acts  to  protect  the  muscle,  (Houk  and 
Hermeman, 1 962) .  The  first  response,  is  Inherently  included  in  the  motor  unit 
recruitment  algorithm  while  the  clasp  knife  reflex  is  modeled  to  reduce  the 
muscle  force  to  zero  should  the  rise  be  greater  than  50  FMAX/sec.  This  value 
is  based  on  maximum  speed  of  motor  unit  recruitment,  ie.  FMAX  in  100  msec 
(Desmedt  and  Godaux,  1977,  1978).  The  limiting  value  of  force  as  t  approaches 
100  msec  (or  n  approaches  1)  is  50  Fp^x/sec. 

D.  Time  Varying  Effects 

The  most  important  time  varying  effect  in  the  muscle  is  fatigue.  It  is 
obvious  that  people  can  maintain  their  maximum  effort  very  briefly  (5 
seconds),  whereas  they  can  maintain  a  force  of  around  a  quarter  of  their 
maximum  strength  for  an  extended  period  of  time.  Such  an  endurance  responses 
can  be  explained  by  examining  the  properties  of  individual  motor  units.  Type 
I  motor  units  tend  to  be  more  aerobic,  less  fatigable  and  are  recruited  at 
lower  tension  levels.  While  Type  II  motor  units  tend  to  be  an  aerobic,  more 
fatigable  and  are  recruited  at  higher  tension  levels,  (Stephens  and  Usherwood, 
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1977).  Although  exact  fatigue  and  recovery  patterns  for  individual  motor 
units  have  not  been  identified,  the  maximum  endurance  time  can  be  estimated 


from  experimental  studies.  The  earliest  experiments  of  Miller  (1932)  implied 
that  the  length  of  time  a  force  could  be  maintained  depended  on  the  fraction 
of  available  strength  to  be  exerted.  This  relationship  was  further  verified 
by  Rohmert,  1960,  Kogi  and  Hakamada,  1962;  Caldwell  1963,  1964;  Monod  and 
Schemer,  1965;  Schutz,  1972.  Only  three  studies  attempted  to  derive  and 
publish  formulas  of  this  relationship.  Monod  and  Schemer  (1965) 
proposed: 
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Kogi  and  Hakamada  (1962)  suggested 


TEND^min^ 


5012 


(*FMAX) 


1.99 


while  Schutz  (1972)  indicated: 
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All  of  these  formulas  have  some  faults  that  limit  their  usefulness  in 
representing  the  empirical  data.  Equation  52  does  not  account  well  for  the 
asymptotic  relationship  of  endurance  approaching  indefinite  times  for  force 
levels  of  15-20JMVC.  Equations  53  and  51*  provide  the  asymptote  but  predict 
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lower  than  normal  endurance  times  for  large  force  levels.  A  separate  formula 
was  developed  for  the  current  work,  based  on  the  data  of  Rohmert  (I960),  who, 
with  over  300  subjects  tested,  had  the  largest  sample  size.  Best  fit  was 
produced  by  the  hyperbolic  relationship  (shown  in  Figure  12): 


TEND(Sec) 


1236.5 


(*W15) 


.618 


-  72.5 


(55) 


Once  the  endurance  time  is  exceeded,  however,  the  person's  strength  does 
not  immediately  fall  to  zero.  For  maximal  or  large  submaximal  efforts,  there 
is  still  a  gradual  decay  to  the  lower  level  of  15~20  percent  found  for 
indefinite  holds  (Petrofsky,  1982,  p.  55).  This  experimental  data  can  be 
modelled  very  easily  using  polynomial  regression  with  time  in  minutes: 

*F  »  98.1  -23. 9t  *  1.9t2  (56) 


V.  PHASE  III  -  MODELLING  THE  GENERAL  MUSCULATURE 

The  attachment  of  the  complete  muscle  systems  to  limb  segments,  includes 
the  identification  of  joint  biomechanics,  the  measurement  of  origin  and 
insertion  coordinates,  the  integration  of  agonist  and  antagonist  actions  and 
computation  of  cross-sectional  areas  for  estimation  of  total  force  production. 
Methods  to  accomplish  this  can  be  best  described  according  to  the  joints  or 
area  of  the  body  involved. 

A,  Elbow  Joint: 

Modelling  of  the  elbow  and  simulation  of  elbow  flexion  is  perhaps  the 
easiest  case  to  examine  and  will  serve  as  a  simple  example  demonstrating  the 
validity  of  the  technique  used  for  the  more  complicated  joints.  The  model 
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Includes  two  segments:  the  humerus,  and  a  combination  of  the  ulna  and  the 
radius  and  three  elbow  flexor  muscles:  biceps  brachii,  brachialis  and 
brachioradialis.  Examination  of  the  biceps  brachii  more  closely  shows  the 
origin  of  the  long  head  to  be  at  or  beyond  the  gleno  -  humeral  joint  (Me  Minn 
and  Hutchings,  1977).  The  insertion  can  be  set  at  approximately  3.5  cm  from 
the  elbow  Joint  corresponding  to  the  data  of  Wilkie  (1950).  Using  the 
cross-sectional  area  of  4.58  cm2  for  the  biceps  (Schumacher  and  Wolff,  1966) 
and  multiplying  by  the  maximum  muscle  force  of  100N/cm2  (Hatze,  1981)  yields  a 
maximum  isometric  tension  of  458  N.  (The  difference  between  anatomic  and 
physiological  areas  due  to  fiber-  orientation  are  accounted  for).  Adding  the 
force  of  690N  generated  by  the  brachialis  (6.9  cm2  x  100  N/cm2)  to  the  biceps 
yields  a  maximum  elbow  flexion  force  of  1148  N.  Such  a  value  can  be  compared 
to  the  data  of  Wilkie  (1950)  who  found  that  his  subjects  could  maintain  a 
maximum  force  of  195.8  N  at  the  wrist.  With  a  lever  ratio  between  muscle 
insertion  distance  and  the  moment  arm  of  the  weight  of  .15,  the  maximum  elbow 
flexion  force  is  1305  N.  These  values  are  remarkably  close  considering  that 
other  minor  muscles  producing  additional  torque  are  not  accounted  for  in  the 
first  calculation.  Similar  calculations  were  conducted  on  the  other  muscles. 

Complete  details  on  all  of  the  elbow  muscles  are  based  or  previously 
collected  data  (Chao  and  Morrey,  1978;  Youm,  et.al.  1979;  Maton,  et.al.  1980; 
Amis,  et.al.,  1980;  Hatze,  1981)  and  are  summarized  in  Table  1,  Appendix  A. 

B.  Shoulder  Joint: 

The  shoulder  is  a  much  more  complicated  joint  consisting  of  three 
separate  joints:  the  glenohumeral  joint,  the  acromioclavicular  joint  and  the 
sternoclavicular  joint.  Correspondingly,  many  more  muscles  are  involved  to 
produce  many  different  actions.  Details  on  the  actions  of  these  muscles, 
points  of  origin  and  insertion,  along  with  cross-sectional  areas  are  taken 
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from  previous  biomechanical  studies  (Schumacher  and  Wolff,  1966;  Dempster, 

1965;  DeLuca  and  Forrest,  1973;  Engin,  1980)  and  given  in  Table  2,  Appendix  A. 

C.  Hip  Joint  and  Knee  Joint: 

Although  consisting  of  only  one  joint,  the  hip  is  a  ball  and  socket  joint 
and  along  with  the  many  muscles  involved,  undergoes  many  different  types  of 
actions.  A  further  complication  is  that  a  majority  of  these  muscles  span  both 
hip  and  knee  joints  and  thus  hip  actions  cannot  be  uniquely  separated  form 
knee  actions.  Details  on  the  actions  of  these  muscles,  points  of  origin  and 
insertion,  along  with  the  cross-sectional  areas  were  taken  from  earlier 
biomechanical  studies  (Merchant,  1965;  Schumacher  and  Wolff,  1966;  Seireg  and 
Arvikar,  1973;  Jensen  and  Davy,  1975;  Crowninshield,  et.al.  1978;  Dostal  and 
Andrews,  1981;  Smidt,  1973;  Nissan,  1980;  Wismanis,  1980;  Minns,  1980)  and  are 
summarized  in  Table  3,  Appendix  A. 

D.  Trunk  and  Neck  Musculature: 

Simulation  of  the  trunk  musculature  is  a  much  more  difficult  undertaking 
than  for  the  previous  joints.  First  of  all,  there  are  many  muscles  involved, 
close  to  20  major  ones  for  the  lower  back  and  trunk  and  equally  many  for  the 
neck  region.  Secondly,  some  of  the  muscles,  such  as  the  longus,  spinalis  and 
semispinalis,  have  many  attachment  sites  between  the  different  vertebrae. 
Thirdly,  the  lines  of  action  of  the  muscle  forces  are  not  always  in  straight 
line,  e.g.,  the  interior  and  exterior  obliques.  Fourthly,  the  vertebral 
Joints  are  complicated  by  the  ligaments  and  their  additional  force-bearing 
capabilities.  Appropriate  approximations  were  used  when  necessary. 

Details  on  insertions  and  origins  were  obtained  from  Rab  et.al.  (1977), 

Rab  (1979),  Takashlma,  et.al.  (1979),  and  Williams  and  Belytschko  (1981),  while 
cross-sectional  areas  were  used  from  Schumacher  and  Wolff  (1966)  and  Williams 
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and  Belytschko  (1981).  A  summary  of  these  findings  is  given  in  Tables  4 
(neck)  and  5( trunk),  Appendix  A. 

VI.  PHASE  IV  -  SOFTWARE  IMPLEMENTATION  OF  ACTIVE  ELEMENTS  INTO  THE  ATB  MODEL 
The  software  modifications  necessary  to  implement  the  active  elements 
into  the  ATB  Model  are  summarized  below.  The  modifications  are  organized  into 
two  sections:  1)  changes  necesary  to  increase  the  number  of  available  muscles 
(belts  within  harnesses),  2)  changes  necessary  to  model  the  active  elements 
within  a  muscle. 

A.  INCREASING  THE  NUMBER  OF  MUSCLES 

As  described  in  a  previous  study  (Freivalds,  1982)  the  incorporation  of 
an  active  neuromusculature  into  the  ATB  model  was  possible  through  the  use  of 
the  advanced  restraint  system  developed  by  Butler  and  Fleck  (1980).  A  harness 
consists  of  one  or  more  belts,  where  each  belt  is  defined  as  a  set  of  straight 
line  segments  connecting  reference  points.  These  points  can  be  placed  within 
the  body  simulating  a  muscle.  Unfortunately  one  of  the  limitations  of  the 
advanced  restraint  system  was  that,  since  being  designed  as  an  external 
harness  system,  a  maximum  of  five  harnesses  could  be  specified.  This  for 
modelling  the  human  musculature  was  an  extreme  limitation  as  was  demonstrated 
by  Freivalds  (1982).  Thus  one  of  the  primary  objectives  was  to  increase  the 
number  of  available  muscles.  This  was  accomplished  by  increasing  the  sizes  of 
several  arrays  and  changing  related  code  necessary  to  handle  these  arrays. 

The  arrays  necessary  for  the  harness/belt  systems  are  found  primarily  in 
common  block  HRNESS.  The  following  arrays  were  increased  in  size  from  20  to 
50: 

XLONG  (50)  -  initial  slacks  for  each  belt 
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NPSTPB  (50)  -  number  of  points  per  belt 
NPTPLY  (50)  -  number  of  points  in  play  per  belt 
NTHRNS  (50)  -  index  to  NTAB  array  defining  the  force  deflection 
functions  for  each  belt. 

Array  NBLTPH  (number  of  belts  per  harness)  was  increased  in  size  from  5  to  50. 
In  common  block  forces  array  BSF  (4,20)  was  changed  to  BSF  (4,50)  to 
accommodate  the  increase  in  the  number  of  belts. 

Using  a  maximum  of  50  muscles  (harnesses)  limits  the  number  of  muscle 
sections  (belts)  to  one  per  harness.  Reducing  the  number  of  muscles 
(harnesses)  still  allows  several  muscles(belts)  per  harness.  Similarly  the 
total  number  of  reference  points  and  points  in  play  for  all  muscles  is  limited 
to  100  (ie.  2  per  each  muscle  for  50  muscles  or  more  if  the  number  of  muscles 
is  reduced). 

In  addition  to  increasing  the  sizes  of  the  arrays  directly  related  to  the 
harness/belt  system,  two  other  arrays  NTAB  and  TAB  in  the  common  block  TABLES 
were  also  expanded.  NTAB  contains  the  index  pointers  to  the  TAB  array  for 
each  function  at  the  allowed  contact  points  and  was  expanded  from  500  to  1250. 
TAB  contains  function  defintions  and  update  information  for  each  contact  point 
and  was  expanded  from  2600  to  6000.  Correspondingly  any  reference  to  variable 
MXNTB,  the  number  of  elements  in  the  NTAB  array,  and  MXNTB2  the  total  number 
of  elements  in  the  TAB  array,  had  to  be  checked.  Consequently  in  subroutine 
FDINIT ,  the  following  lines  were  altered: 

Line  #  51  IF  (MXTB2.GT.5000)  WRITE  (6.62)  MXTB2) 

Line  #54  IF' (MXNTB. GT. 1250)  WRITE  (6.63)  MXNTB 

Line  #52  IF  (MXTB2.GT.5000  OR  MXNTB. GT . 1 250 )  STOP  16 

Cases  of  the  above  arrays  being  initialized  or  set  equivalent  to 


other  arrays  were  also  checked.  Thus  in  subroutine  OUTPUT  in  order  to 
initialize  the  arrays  in  common  block  FORCES,  the  limit  on  the  counter  in  the 
DO  loop  in  line  #44  must  be  increased  from  760  to  880: 

Line  #44  DO  1 1  1-1,880 

In  subroutine  RSTART,  new  values  of  variables  to  be  input  are  set 
equivalent  to  labelled  common  blocks.  Thus  the  following  lines  were  altered. 

Line  #76  DIMENSION  RC1 1  (880),  IC1 1  (9) 

Line  #93  DIMENSION  IC14  (1304) 

Line  #139  DIMENSION  RC23  (1952),  IC23  (1050) 

In  subroutine  SEARCH  called  by  subroutine  RSTART  to  compute  the  labelled 
common  block  numbers  and  items,  corresponding  changes  again  had  to  be  made. 
Thus  the  following  lines  were  altered  to  become: 

Line  #7  DIMENSION  BVAR  (260),K0UNT  (25),  NDIM  (3,260),  NJ(3),  NK(3), 

INDEX( 3 ) 

Line  #30  DIMENSION  C2305),  NC23(45) 

Line  #55  EQUIVALENCE  (C23(1),  BVAR(241)),  (NC23(1),  NDIM( 1 ,241 ) ) . 

Line  #56  EQUIVALENCE  (C24(1),  BVAR  (253)),  (NC24(1),  NDIM0.253)) 

Line  #58  DATA  NVAR/260/.K0M/24/,  BLANK/8H  / 

Line  #152  DATA  NC1 1  /7,300,  4,  50,  0,  10,  20,  0,  3,  20,  0,  7,  30,  0, 

Line  #176  1250,  0,  0,  6000,  0,  0/ 

Line  #240  add,  8H  FSCALE  ,  8H  FTYP1  ,  8H  FMU  / 

Line  #242  2,  0,  0,  5,  100,  0,  50,  0,  0,  50,  0,  0, 

Line  #243  50,  0,  0,  50,  0,  0,  50,  0,  0,  50,  0,  0,  50,  0,  0/ 

To  correct  for  the  additional  output  arising  from  a  large  number  of 
muscle/harnesses,  one  additional  modification  may  be  needed.  Should  the  user 
desire  tabular  time  histories  of  the  various  events.  NPRT  (4)  on  input  an  A5. 
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must  be  set  according  to  the  specifications  in  the  User's  Manual  (Fleck  and 
Butler,  1982).  The  Tabular  time  histories  can  be  of  two  types:  1)  optional 
output  controlled  by  input  Cards  H.1  to  H,7  and  2)  automatic  output.  The 
optional  tabular  time  histories  are  as  follows: 

la)  linear  accelerations  for  specified  points 

lb)  linear  velocities  for  specified  points 

lc)  linear  displacements  for  specified  points 

ld)  angular  accelerations  for  specified  points 

le)  angular  velocities  for  specified  points 

lf)  angular  rotations  for  specified  points 

lg)  joint  torques 

The  automatic ‘tabular  time  histories  are  as  follows: 

2a)  plane-segment  contacts 
2b)  belt-segment  contacts 
2c)  muscle-harness  forces 
2d)  spring  damper  forces 
2e)  segment-segment  contacts 
2f )  airbag-segment  contacts 

During  program  execution  the  above  calculated  values  are  stored  in 
the  array  ZTTH  in  common  block  TMPVS  by  subroutine  POSTPR  and  printed  by 
subroutine  HEDING.  In  the  present  version,  array  ZTTH  was  dimensioned 
(1*J, **5, 2)  allowing  theoretically  approximately  7  tabular  time  histories. 
However,  since  the  common  block  TEMPVS  was  dimensioned  to  a  larger  size  in  the 
main  program  array,  ZTTH  overran  its  normal  boundaries  but  stayed  within  the 
size  of  the  total  common  block  TEMPVS  which  is  equal  to  10538  integer  words 
By  default  this  allowed  approximately  38  tabular  time  histories.  With  an 


increased  number  of  muscle/harness  systems  and  corresponding  automatic  tabular 
time  histories  there  is  much  greater  chance  to  overflow  common  block  TEMPUS 
and  crash  the  program.  Thus  array  JTMPVS  in  common  block  TEMPVS  in  the  main 
program  was  increased  from  10538  to  40,000  integer  words.  (Each  integer  word 
is  equal  to  four  bytes).  This  should  allow  approximately  200  tabular  time 
histories  which  should  be  sufficient  for  50  muscle/ harness  systems  and  a 
variety  of  other  segment-plane-airbag  interactions.  In  case  there  is  a 
problem  with  array  size  limitations  and  overflow,  the  dimension  of  array 
JTMPVS  can  be  increased  according  to  the  number  of  points  or  contacts 
specified  in  1a)-lg)  and  2a)-2f)  above.  Thus 
//  OPTIONAL  ( #1 a+//1 b+#1 c+Z/1 d+1 e+1 f )  X  4  +  (#1g)  X  7 

#  AUTOMATIC  -  (#2a)  X  7  +  (#2 b  +  #2c  +  #2d)  X  4  +  (#2e)  xIO  +  (#2f)  x  3 
//  INTEGER  WORDS  FOR  JTMPVS  =  (//OPTIONAL  +  //AUTOMATIC)  X  45  +  3560 

B.  ACTIVE  NEUROMUSCULATURE 

As  described  previously  in  earlier  works  (Freivalds,  1982;  Freivalds, 
1984),  the  incorporation  of  active  muscle  elements  was  accomplished  through 
the  use  of  the  advanced  restraint  system  developed  by  Butler  and  Fleck 
(1980)  . 

In  the  present  effort,  instead  of  merely  manipulating  input  data  and 
using  the  existing  software  to  simulate  an  active  neuromusculature,  an 
additional  subroutine  named  MUSCLE,  was  written  and  incorporated  into  the  ATB 
Model.  Subroutine  MUSCLE,  however,  still  maintains  the  basic  or 
organizational  relationships  of  the  active  restraint  system,  but  only 
specifies  more  accurately  active  muscle  parameters.  This  modification  is 
transparent  to  the  user  allowing  the  flexibility  of  using  old  input  data 
without  musculature  as  well  as  new  input  data  with  musculature  on  the  same 
software. 
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The  modifications  needed  to  incorporate  an  active  neuromusculature  in  the 
present  ATB  Model  can  be  summarized  as  follows: 

1)  A  new  subroutine  named  MUSCLE  was  written  and  incorporated  into  the 
code 

2)  Major  changes  occurred  in  subroutine  HINPUT 

3)  Minor  changes  occurred  in  subroutines  FDINIT  and  HBELT. 

These  modifications  will  be  described  in  detail  in  the  following 
sections . 

1.  Subroutine  MUSCLE 

The  complete  code  for  subroutine  MUSCLE  is  given  in  Appendix  B  and  a 
flowchart  of  the  logic  is  given  in  Fig.  13.  The  basic  logic  is  as  follows. 
Upon  entry  to  subroutine  MUSCLE  the  current  time  is  adjusted  by  the  delay  time 
which  is  specified  by  the  user  in  subroutine  HINPUT.  Thus  the  user  can  insert 
a  10-30  msec  delay  to  account  for  nerve  conduction  times,  a  150-200  msec  delay 
to  account  for  reaction  times,  or  even  longer  delays  depending  on  the 
circumstances.  Irregardless  of  the  delay  time,  the  passive  (spring  and 
viscous  damping)  muscle  forces  are  calculated  continuously.  The  active  muscle 
forces  are  calculated  only  when  the  current  time  exceeds  the  delay  time. 

Next  the  percent  of  motor  units  ($  mu)  recruited  for  that  given  time  is 
calculated.  The  rate  of  recruitment  is  specified  by  the  user  in  subrouting 
HINPUT.  This  rate  can  vary  from  a  maximum  rate  of  100$  mu  recruited  (or 
maximum  force)  in  100  msec  to  much  slower  rates  in  seconds.  The  program 
checks  that  this  maximum  rate,  based  on  the  ballistic  muscle  contraction  data 
collected  by  Desmedt  and  Godaux  (1977,  1978),  is  not  exceeded. 

Next  the  program  compares  the  $  mu  recruited  with  the  percent  of  Type  I 
mu  as  specified  by  the  user  (in  subroutine  HINPUT).  If  the  %  mu  recruited  is 
less  than  the  %  of  Type  I  mu  then  contraction  times  based  on  Type  I  mu 


40 


properties  are  calculated.  If  the  %  mu  recruited  exceeds  the  %  of  Type  I  mu, 
then  contraction  times  based  on  Type  II  mu  properties  are  used.  Next  the 
fractional  forces  developed  by  the  various  elements  as  described  in  Chapter  I 
are  calculated.  The  total  muscle  force  based  on  %  mu  recruited  is  then 
calculated  and  adjusted  by  the  fractions  contributed  by  each  element.  Finally 
an  endurance  time  based  on  the  total  muscle  force  is  calculated.  If  the 
current  time  is  less  than  the  endurance  time,  the  muscle  force  remains 
unaltered.  However,  if  the  current  time  is  greater  than  the  endurance  time 
the  muscle  force  is  reduced  proportionately  according  to  Eq.  56. 

2.  Subroutine  HINPUT 

The  complete  code  for  subroutine  HINPUT  is  given  in  Appendix  C.  The 
major  changes  are  involved  in  supplying  additional  muscle  parameters  on  the 
F.8.C  card.  These  changes  occur  between  lines  53-57  in  subroutine  HINPUT. 

The  primary  input  in  line  # 53  becomes, 

READ  (5,  14)  NF,  XLONG  (J) ,FSCALE( J)  FTYPI(J),  FMU(J),  RECTIM(J),  DEL(J). 

Where  FSCALE(J)  »  Scaling  factor  for  maximum  muscle  force 

FTYPI(J)  -  percent  Type  I  motor  units  (range:  0-100) 

FMU(J)  »  percent  motor  units  recruited  (range:  0-100) 

RECTIM(J)  -  Motor  unit  recruitment  time  (seconds) 

DEL(J)  -  delay  time  before  motor  unit  recruitment  (seconds) 
and  the  format  in  line  #54  becomes: 

14  FORMAT  (514,  F12.0,  5F8.0) 

Thus  the  format  and  order  for  the  previous  input  variables  remains  the  same 
and  the  additions  become  transparent  to  a  user  not  concerned  with  them.  Since 
the  E  card  functions,  previously  used  to  define  the  muscle  characteristics, 
are  not  needed,  these  can  be  left  blank  or  zeroed.  However  NF(1)  must  be 
specified  with  a  negative  value  to  indicate  that  a  muscle  and  not  a  harness  is 
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being  defined.  Since  the  previous  version  allowed  only  NF(2),  NF(3)  and  NF(4) 
negative  for  rate  dependent  functions,  but  not  NF(1),  these  become  mutually 
exclusive  definitions. 

Then  depending  on  the  value  of  NF(1),  two  different  echos  are  produced. 

If  NF(1)  £0  the*  previous  format  in  line  #55  is  left  unaltered: 

IF  (NF(I).GE.O)  WRITE  (6,15)  I,  J,NF,XL0NG(J)  .(JNITL 
whereas  for  the  current  version  and  NF(1)  S  0  the  following  is  inserted: 

IF  (NF(I).LT.O)  WRITE  (6,19)  I,  J.FSCALE  (J),  UNITM,  FTYPI(J),  FMU(J), 
RECTIM(J) ,  DEL(J) 

with  its  aproprlate  format  statement.  Also,  several  statements  are  inserted 

« 

to  check  whether  the  input  parameters  are  within  limits. 

1)  The  percent  of  Type  1  mu  must  be  less  than  100: 

IF  (FTYPI(J).GE.  100.)  FTYPI(J)  -  99.99 

2)  The  percent  of  Type  1  mu  must  be  greater  than  or  equal  to  zero. 

IF  (FTYPI(J).LE.O)  FTYPHJ)  -  .0001 

3)  The  first  30%  of  mu  recruited  must  be  Type  I. 

IF  (FMU(J) . LE. 30.  .AND.  FMU( J) .GT. FTYP1 ( J) )  FTYPI(J)  -  FMU(J) 

4)  Finally  two  statements  change  percent  values  to  fractional  values. 

FTYP1 ( J)  -  FTYP1 ( J)/100. 

FMU(J)  -  FMU ( J) /1 00 . 

Revised  instruactions  for  defining  the  F.8.  cards  are  given  in  Appendix  D. 

3.  Minor  Modifications 

In  subroutine  HBELT  at  line#84,  depending  on  the  value  of  the  first 
harness  function,  ie.  NF(1),  either  the  new  subroutine  MUSCLE  or  the  old 
subroutine  FRCDFL  is  called  to  calculate  the  forces  due  to  strain  and  strain 
rate  thus: 

NM  -  NTA8  ( NT  + 1 ) 
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IF  (NM.LT.O)  CALL  MUSCLE  (etc) 

IF  (NM.GE.O)  CALL  FRCDFL  (etc) 

Also  Line  #78  was  eliminated  to  allow  for  negative  stains,  which  is  necessary 
for  concentric  muscle  motion. 

In  subroutine  FDINIT,  to  eliminate  the  need  of  storing  E-card  functions 
not  required  for  the  active  musculature,  three  lines  were  inserted  after  line 
#21  : 

IF  ( NF ( 1 ) . GT . 0 )  GO  TO  50 
NTAB(NT)  -  -1 
GO  TO  56 

and  line  #22  was  given  line  number  50. 

VII.  SIMULATION  AND  VALIDATION 

A.  SIMULATIONS  OF  TRUNK  AND  NECK  MUSCULATURE 

In  Phase  IV,  the  advanced  neuromuscular  model  was  validated  via  simulation 
of  human  body  responses  to  high  G( lateral)  acceleration.  Data  obtained  under 
similar  conditions  on  air  crew  personnel  experiencing  tests  in  the  Dynamic 
Enviornmental  Simualtor  at  Wright-Patterson  AFB,  OH,  were  used  for  comparison 
purposes.  Some  of  this  data  was  used  in  an  earlier  pilot  study  developing  a 
trunk  musculature  (Freivalds,  1984). 

The  full  ATB  Model  with  15  body  segments  (head,  neck,  upper  torso,  center 
torso,  pelvis,  upper  arms,  lower  arms,  upper  legs,  lower  legs, and  feet)  and  14 
joints  (head  junction,  neck  junction,  waist,  L5/S1 ,  joint,  two  hips,  two  knees 
two  ankles,  two  shoulders,  two  elbows)  was  utilized  to  provide  an  adequate  human 
neuromusculature.  The  muscles  described  in  Chapter  V  and  summarized  in  Tables 
1-5  needed  to  be  added  to  the  existing  ATB  Model.  For  each  muscle  specific 
coordinates  of  origin  and  insertion  points  were  determined  form  anatomical  texts 
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(Qulrlng  et.al.,  1945;  Quiring,  1 947 ;  Gray,  1974;  McMinn  and  Hutchings,  1977). 
Cross  sectional  areas  obtained  from  Schumacher  and  Wolff  (1966)  and  Williams  and 
Belytschlso  (1981)  were  multiplied  by  the  muscular  force  constant  of  100  N/cm^ 
(Hatze,198l)  to  determine  a  muscle  force  scaling  factor.  These  resulting  values 
were  then  converted  to  English  units  for  use  in  the  present  ATB  Model  mode  and 
are  summarized  in  Tables  6-10,  Appendix  E. 

For  the  simulations  of  responses  to  lateral  G  forces,  the  body  segments 
were  arranged  in  the  semi-reclining  posture  maintained  by  air  crew  personnel  in 
the  cockpit.  The  lower  trunk  was  restrained  by  a  lap  belt;  the  feet  fixed  on 
the  floor;  any  other  restraints,  such  as  shoulder  pads  or  hands  placed  on 
controls,  were  eliminated.  Only  the  neck  and  trunk  musculature  on  the  right 
side  of  the  body  were  activated  so  as  to  reduce  program  complexity  and  execution 
times.  A  2  Gy  lateral  force  was  applied  to  the  body  and  the  acceleration; 
velocity  and  displacement  of  various  body  segments  were  recorded. 

A  graphical  response  of  the  whole  body  response  (using  only  trunk  and  neck 
musculature)  to  the  lateral  force  over  time  is  shown  in  Figure  14.  For 
comparison  purposes,  the  response  to  a  control  case  with  no  musculature  is  given 
in  Figure  15a,  while  the  response  using  the  previous  simplified  musculature  is 
given  in  Figure  15b.  Although  the  musculature  does  not  completely  prevent  the 
lateral  deflection  of  the  body,  the  response  is  significantly  delayed  with  head 
and  neck  maintaining  the  upright  position  for  a  longer  period  of  time.  The 
result  is  better  observed  in  Figure  16,  which  shows  the  plot  of  angular 
displacement  of  the  upper  trunk  for  all  three  conditions.  At  the  end  of  256 
msec  the  angular  displacement  is  reduced  by  20°  with  the  use  of  musculature. 
Based  on  the  time  history,  the  response  with  the  musculature  lags  up  to  *40  msec 
behind  the  control  response 

B.  SENSITIVITY  OF  MUSCLE  ORIGINS  AND  INSERTIONS 
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Control  condition 


Graphical  Display  of  Trunk  Response  to  Lateral  G  Forces  (2G) 


ANGULAR  DISPLACEMENT  (ROLL)  OF  UPPER 
TRUNK  DUE  TO  LATERAL  G-FORCES 


TIME  -  msec 


One  Important  factor  effecting  the  accuracy  in  stimulating  the  active 
neuromusculature  is  the  accuracy  of  the  placement  of  the  musculature,  i.e.  the 
origins  and  insertions  of  the  muscles  being  used.  The  origins  and  insertions 
specified  in  Appendix  E  were  derived  from  anatomy  texts  such  as  McMinn  and 
Hutchings  (1977),  typically  by  measuring  directly  from  scale  photographs  of  the 
anatomical  structures.  Obviously,  some  human  error  is  involved,  especially  in 
using  two  dimensional  photographs  of  three  dimensional  structures.  Thus  it  was 
decided  to  simulate  errors  in  muscle  placement  and  compare  the  resulting  muscle 
forces  with  "true”  muscle  placements. 

Five  different  cases  were  used  in  simulating  elbow  flexion  with  three 
muscles:  the  biceps  brachii,  the  brachialis  and  the  brachloradialis.  In  Case 
1,  the  origin  of  each  muscle  was  displaced  by  1 0%  of  the  insertion  distance.  In 
Case  3  both  origins  and  Insertions  were  shifted  by  1 0% .  In  Case  4  both  origins 
and  insertions  were  shifted  by  20$  of  the  distance.  Case  5  was  the  control 
condition  with  "true"  origins  and  insertions.  The  simulation  was  of  elbow 
extension  using  a  weight  of  60  lbs.  held  in  the  hand  as  shown  in  Fig.  17.  Since 
this  weight  exceeded  the  maximum  weight  that  could  be  maintained  at  90°  of  elbow 
flexion,  the  muscle  were  were  forcibly  extended  in  an  eccentric  fashion.  The 
resulting  muscle  forces  were  plotted  as  a  function  of  the  Included  elbow  angle 
in  Figs.  18-20. 

The  general  motion  of  the  curves  can  be  best  explained  by  the  forces 
produced  by  each  of  the  elements  Involved  in  the  muscle  model.  The  rather  sharp 
overall  increase  is  due  to  the  force  velocity  relationship  for  eccentric  motion. 
The  damping  function  also  adds  some  force  with  increasing  speed.  But  as  the  arm 
is  extended  further,  the  length  tension  relationship  contributes  significantly 
in  reducing  the  force.  The  brachioradialis  is  most  affected,  the  biceps  the 
least.  However,  for  all  four  cases  there  is  minimal  deviation  from  the  norm 
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Fig.  17  ANATOMY  OF  THE  ELBOW  FLEXORS 


FORCE  -  lbs. 


BICEPS  BRACHII 


INCLUDED  ANGLE  -  degrees 


Fig.  20 


Force  Production  by  the  Biceps  during  Eccentric  Hoi  ion 


with  an  average  deviation  of  about  3 %.  The  maximum  deviation  is  9$.  Thus  for 
most  practical  purposes,  deviations  in  muscle  origins  and  insertions  will  not  be 
significant  for  eccentric  motions. 

Concentric  muscle  contractions,  using  a  weight  of  15  lbs.,  were  also 
simulated.  The  results  are  plotted  an  described  in  detail  in  Appendix  F.  In 
general,  for  all  the  curves,  there  was  a  sharp  decrease  in  force  due  to  the 
force  velocity  relationships,  with  a  later  recovery  as  the  velocity  decreased. 
There  were  greater  deviations  from  the  norm  -  averging  about  5t  but  reaching  up 
to  10$.  Again  the  brachioradiales  showed  greatest  effect  especially  for  Case  2. 
This  is  due  to  a  configuration  opposite  of  the  biceps  and  brachialis,  acting 
more  as  a  shunt  rather  than  a  spurt  muscle  (MacConail,  1949). 

Only  a  slight  deviation  is  noticable  for  Cases  3  and  4.  Perhaps  this  can 
be  rationalized  by  having  the  effects  of  both  offsets  cancel  each  other  out. 

C.  ANATOMICAL  STUDIES  OF  TRUNK  MUSCULATURE 

Other  simulations  of  trunk  musculature  were  also  attempted.  Special 
efforts  were  applied  to  generate  accurate  crossectlonal  areas,  as  well  as  muscle 
origin  and  insertion  points.  The  details  of  muscle  cross  sectional  areas  as 
Investigated  by  C.  S.  Davis  are  presented  in  Appendix  A. 

vm.  CONCLUSIONS 

Further  simulations  of  various  human  body  motions  or  responses  to  external 
forces  need  to  be  conducted  in  order  to  more  adequately  validate  the  model. 

These  would  include  not  only  additional  lateral  forces,  but  also  forward/back 
(Gx)  and  up/down  (Gz)  forces  as  well  as  multidirectional  forces.  Various 
combinations  of  muscle  parameter  values  need  to  be  tested,  so  as  to  determine 
the  optimum  values  for  generating  the  most  realistic  human  response. 
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Similarly,  further  examination  of  other  neuromuscular  reflexes  may  still  lead 
to  further  improvements  in  the  model.  However,  the  simulations  demonstrated 
that  the  presently  developed  muscle  model  can  adequately  represent  an  active 
human  neuromusculature  response  to  dynamic  mechanical  stresses  and  can  serve 
as  a  cost  effective  research  and  developmental  tool. 
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APPENDIX  A 


INPUT  DATA  FOR  HUMAN  MUSCULATURE 


Table  1.  Specifications  on  Elbow  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2! 

Biceps  brachii 

Short  head  from 
coracoid  process 
of  scapula,  long 
head  from 
supraglenoid 
tuberosity 

Radial  tuberosity 

Flexion  of 
forearm 

3.55 

Brachialis 

Lower  anterior 
surface  of 
humerus 

Coronoid  tuber¬ 
osity  of  ulna 

Flexion  of 
forearm 

4.63 

Brachiora- 

dialis 

Proximal  two- 
thirds  of 
humerus 

Styloid  process 
of  ulna 

Flexion  of 
forearm 

1.37 

Triceps 

brachii 

Long  head  from 
infraglenoid 
tunerosity  of 
scapula,  lateral 
and  medial  head 
from  posturior 
surface  of  humerus 

Olecranon 

Extension  of 
f  orearm 

16.38 

Anconeus 

Lateral  epicondyle 
humerous 

Olecranon 

Extension  of 
forearm 

.94 

Pronator 

tere3 

Medial  epicondyle 
of  humerus 

Middle  of  radius 

Pronation 
of  forearm 

1 .61 

Supinator 

Lateral 
epicondyle  of 
humerus 

Laterial  and 
anterior  surface 
of  radius 

Supination 
of  forearm 

1.77 

Table  2.  Specifications  on  Shoulder  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2) 

Deltoid 

Clavicle,  scapula 
acromion 

Deltoid  tuberosity 
of  humerus 

Abduction 
of  arm 

!  11.01 

Supraspinatus 

Supraspinous  fossa 
of  scapula 

Greater  tuberosity 
of  humerus 

Abduction 
of  arm 

3.3 

Pectoralls 

major 

Clavicle,  Stermum 

Bicipital  groove 
of  humerus 

Adduction 
of  arm 

6.8 

Latisslmus 

dorsi 

Lower  Thoracic 
and  lumbar 
vertebrae 

Bicipital  groove 
of  humerus 

Adduction 
of  arm 

5.37 

Teres  major 

Inferior  angle  of 
scapula 

Bicipital  groove 
of  humerus 

Adduction 
of  arm 

4.97 

Teres  minor 

Axillary  border 
scapula 

Greater  tuberosity 
of  humerus 

Adduction 
of  arm 

1 .57 

Subscapularis 

Subscapular  fossa 
scapula 

Lesser  tuberosity 
of  humerus  j 

Flexion 
of  arm 

9.9 

Coracobrachi 

alls 

Coracoid  process 
scapula 

Medial  border  of 
humerus  ! 

1 

Flexion  of 
arm  ! 

1.52 

Infraspinatus 

Infrasplnous 
fossa  of  scapula 

Greater  tuberosity 
of  humerus 

Extension 
of  arm 

5.98 

Table  3  Specifications  on  Hip  and  Knee  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion  i 

Action 

Area  (cm^) 

Gluteus 

medius 

Iliac  crest  of 
pelvis 

Greater  trochanter 
of  femur 

Abduction 
of  thigh 

21  .18 

Gluteus 

minimus 

Outer  surface  of 
ilium  (pelvis) 

Greater  trochanter 
of  femur 

Abduction 
of  thigh 

9.6 

Tensor  fasiae 
latae 

Anterior  part  of 
iliac  crest  of 
pelvis 

Iliotibial  tract 
of  femur 

Abduction 
flexion  of 
thigh 

2.48 

Obturatus 

internus 

Obturatus  foramen 
area  of  pelvis 

Greater  trochanter 
of  femur 

Abduction 
of  thigh 

3.91 

Adductor 

longus 

Pubis  of  pelvis 

.  _J 

Linea  aspera  of 
femur 

Adduction 
flexion  of 
thigh 

5.03 

Adductor 

brevis 

1 

Pubis  of  pelvis 

Below  lesser 
trochanter  of 
femur  ! 

Adduction 
flexion  of 
thigh 

4.54 

Adductor 

magnus 

Ischial 
tuberosity  of 
pelvis 

Linea  aspera  of 
femur 

Adduction 
flexion  of 
thigh 

20.58 

Pecti neus 

Pubic  tubercle 
of  pelvis 

Below  lesser 
trochanter  of 
femur 

Adduction 
flexion  of 
thigh 

2.47 

Quadratus 
f  emor i s 

Ischial 
tuberosity  of 
pelvis 

Quadrate  tubercle 
of  femur 

Adduction 
of  thigh 

2.91 

Obturatus 

externus 

Obturator  foramen 
area  of  pelvis 

Trochanteric 
fossa  of  femur 

Adduction 
of  thigh 

4.95 

Gluteus 

maximus 

Iliac  crest  of 

sacrum 

Iliotibial  tract 
of  femur 

Extension 
of  thigh 

29.42 

Semimem¬ 

branosus 

Ischial  tuberosity 
of  pelvis 

Upper  part  of 
tibia 

Extension 
of  thigh 
flexion  of 
leg 

12.97 

Semitendi- 

nosus 

Ischial  tuberosity 
of  pelvis 

Medial  condyle 
tibia 

Extension 
of  thigh 

4.33 

Table  3  (Continued) 


Biceps 
f emoris 

Ischial 
tuberosity  of 
pelvis,  linea 
aspera  of  femur 

Lateral  condyle 
of  tibia,  head 
of  fibula 

Extension 
of  thigh, 
flexion  of 
leg 

Quadriceps 
f emoris 

Iliac  spine  of 
pelvis  anterion 
surface  of 
femur 

Patela 

Flexion  of 
thigh 
extension 
of  leg 

56.0 

Iliopsoas 

L2-L4  vertebral 
bodies,  iliac 
fossa  of  pelvis 

Lesser  trochanter 
of  femur 

Flexion  of 
thigh 

15.06 

Gastrocnemus 

Medial  and 
lateral  condyles 
of  femur 

Calcaneus 

Flexion  of 
leg 

15.66 

Popliteus 

Lateral  condyle 
of  femur 

Posterior  surface 
of  tibia 

Flexion 

and 

rotation 
of  leg 

1.99 

Gracilis 

! 

Pubic  symphysis 
of  pelvis 

Upper  medial 
surface  of 
tibia 

Flexion  of 
leg, 

adduction 
of  thigh 

1 .63 

Sartor ius 

Iliac  Notch  of 
of  pelvis 

Upper  Medial 
surface  of 
tibia 

Flexion  of 
leg  and 
thigh 

1 .55 

Table  4  Specifications  on  Neck  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2) 

Longus  Capitis 

Transverse 
processes  of  C3, 
C4 ,  C5,  C6 

Basilar  part  of 
occipital  bone 

Flexion  of 
head 

.75 

Rectus  capitis 
interior 

Cl  transverse 
processes 

Front  of  foramen 
magmum  on 
occipital  bone 

Flexion  and 
notation  of 
head 

.25 

Rectus  capitis 
lateralis 

Cl  transverse 
processes 

Jugular  process 
of  occipital 
bone 

Lateral 
flexion  of 
head 

.25 

Rectus  capitis 
posterior  major 

C2  spinous 
process 

Inferior  nuchal 
line  of  occipital 
bone 

Extension, 
lateral 
flexion  of 
head 

.50 

Rectus  capitis 
posterior  minor 

Cl  spinous 
process 

Inferior  nuchal 
line  of  occipital 
bone 

Extension, 
laterial 
flexion  of 
head 

.385 

Obliquus 

capitis 

superior 

Cl  transverse 
processes 

Inferior  nuchal 
line  of  occipital 
bone 

Extension 
and  lateral 
rotation  of 
head 

1.00 

Splenius 

capitis 

Spinous  process 
of  T1  and  C7 

Occipital  bone 
and  temporal 
bone 

Extension 
and  lateral 
flexion  of 
head 

1  .22 

Longis3imus 

capitis 

Transverse 
processes  of 

T1,  C6,  C4 

Mastoid  process 
of  Temporal 
bone 

Extension 
and  Lateral 
flexion  of 
head 

.5 

Spinalis 

capitis 

Transverse  pro¬ 
cesses  of  T1  and 
C7 

Between  superior 
and  Inferior 
nuchal  line  of 
occipital  bone 

Extension 
and  lateral 
flexion  of 
head 

.5 

Semispinalis 

capitis 

Transverse  pro¬ 
cesses  of  T1  and 
C7 

Between  superior 
and  inferior 
nuchal  line  of 
occipital  bone 

Extension 
and  lateral 
flexion  of 
head 

2.38 
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Table  M  (Continued) 


Trapezius 


Levator  scapula 


Scalenus 
anterior,  mediu 
and  posterior 


Splenius 

cervicis 


Longissimus 

cervicis 


Spinalis 

cervicis 


Semispinalls 

cervicis 


Heads  of 
clavicles  and 
and  spines  of 
scapulae 

Occipital  bone 
and  thoracie 
vertebrae 

-  Head  of  sternum, 
medial  sections 
and  heads  of 
clavicles 

Occipital  bone 
and  temporal 
bone 

e  Medical  sections 
of  scapulae 

Transverse  pro¬ 
cesses  of  Cl ,  C3 

Anterior  side  of 
body  of  C5 
anterior  side  of 
body  of  C6 
anterior  side  of 
body  of  T1 

Anterior  side  of 
body  of  CM 
anterior  side  of 
body  of  C3 
anterior  side  of 
body  of  CM 

Medial  clavicle 

s 

C3,  CM,  C5.  C6, 

C7  transverse 
processes 

T1  spinous  pro¬ 
cess 

Transverse  pro¬ 
cesses  of  Cl ,  C2 

C3 

T1  transverse 
processes 

Transverse  pro¬ 
cesses  of  C2,  C3, 
CM ,  C5,  C6 

Spinous  process 
of  T1  and  C7 

C2  spinous  process 

Transverse  pro¬ 
cesses  of  T2  and 
C7 

Spinous  process 
of  C2,  C3,  CM  and 
C5 

Extension 
lateral 
flexion  of 
head 


Lateral 
flexion  and 
flexion  of 
head 


Lateral 
flexion  of 
neck 


Flexion  of 
neck 


Flexion  and 
laterlal 
flexion  of 
neck 


Extension 


Table  5  Specifications  on  Trunk  Musculature 


Cross- 

Sectional 


Muscle  Group 

Origin 

Insertion 

Action 

Area  (cm2) 

Iliocastalis 

lumborum 

Transverse  pro¬ 
cesses  L1-L5 

- 

Lower  six  ribs 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1 .0 

Iliocostalis 

dorsi 

Lower  six  ribs 

Upper  six  ribs 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

.5 

Longissimus 

dorsi 

Transverse  pro¬ 
cesses  L1-L5 

Transverse  pro¬ 
cesses  T1-T12 

1 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1 .0 

Spinalis  dorsi 

Spinous  processes 
L2,  LI,  T12,  Til 

( 

Spinous  processes 
T4-T8 

. 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1  .0 

Semispinalis 

dorsi 

Transverse  pro¬ 
cesses  T7-T12 

Spinous  processes 
C6,  C7,  T1 ,  T2, 

T3 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

1 .0 

Multif idus 

i 

Transverse 
processes  C5-T12 

| 

Spinous  processes 
above  vertebra  of 
origin 

Extension 
and  lateral 
flexion  of 
vertebral 
column 

Interspinales 

Spinous  pro¬ 
cesses  L5-C2 

Spinous  processes 
above  vertebra 
of  origin 

BM 

.5 

Intertrans- 

versaril 

Connect  adjacent 

transverse  processes 

Lateral 
flexion  of 
vertebral 
column 

.25 

Table  5  (Continued) 


External 

Oblique 

Anterior  half  of 
iliac  crest 

Lower  8  ribs 

Compresses 

abdomen 

6 . 85 

Internal  oblique 

■  Anterior  half  of 
iliac  crest 

Lower  3  ribs 
and  mid  line  of 
body 

Compresses 

abdomen 

5.68 

Rectus 

abdominus 

Pubic  symphysis 

Xyphoid  process 

Flexes 

vertebral 

column 

2.66 

Quadratus 

Lumborum 


Iliac  crest 


Transverse  pro¬ 
cesses  L1-L4 


Flexes 

vertebral 

column 
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APPENDIX  B 

LISTING  OF  FORTRAN  IV  SOURCE  DECK  FOR 
SUBROUTINE  MUSCLE 


SUBROUTINE  MUSCLE  ( 0, RATE, M,N,F, GLOSS. IM.XTIME) 

C  REV  6/30/82 

IMPLICIT  REAL*8( A-H.O-Z ) 

CO'HON/TAHLES/MXNTI  ,MXNTli,MXTBI  ,'<xrB2,Mri  (BO),NTAB(  1  2‘jO)  ,  TAB(  6000) 
COM.WON/HMNESS/  BAR(  1 2 , 1  00) ,  BB(1  DO) ,  BUDOT ( I  00) , PLOSS  (2 ,  I  00), 

*  XL0N0(S0),HTIMr:(2) ,  I  BARIS,  1 00)  ,NL(2,  I  OJ  ) . 

*  NPTSPU(bO) , MpTPLY (jO  )  ,NTHRNS(bO>  JJoLTPHtoO) 

*  ,  FSCALEf  50 ) ,  FTYP I  ( 20  )  ,FMU  (SO  )  ,,:ECT  I M  (oo  )  ,[)EL<  oU) 
DIMENSION  FPREV(bO) .TPREV(bO) 

F  =  0.0 
ELOSS  <=0.0 
IFfN.NE.  I  )G0r099 
C  CALCULATE  MUSCLE  FORCES 
C 

C  ACCOUNT  FOR  REACTION  OR  OTHER  TIME  DELAY 
TIME=xriME-DEL( I M ) 

IFCTIME.LE.O. )TIME=0. 

C  MU  RECRUITED  LINEARLY  FOR  FIRST  RECTIM  SEC 
FN-TI ME/RECTI M( Im) 

IF<  FN.LE.O.O)FN».OOl 

C  CHECK  IF  f.,U  RECRUITMENT  LEVEL  IS  EXCEEDED 
I  Ft  FN.GE.FMUdM  ))  FN»FMU<  IM) 

C  CHECK  IF  MU  RECRUITMENT  EXCEEDS  FRACTION  OF  TYPE  I  MU 
I F( FN.GE.FTYpI ( IM )) GOTO 60 
C  CALCULATE  AVERAGE  CONTRACTION  TIMES 
E5=EXP(5.*FN)*(b.*FN-l .  ) 

TC=. I-.OI l*(ED+l .)/FTYPI(IM)/(EXP(2.*FN)-l . ) 

Goro7o 

60  EbF=EXP(b.*r'TYPI  ( IM  ) ) *( b. *FTYP I  ( IM)-I  . ) 

TC  l*>.  i  -  .01 1  ★(  EbF+ 1  .  )/FTYPl(IM)/(EXP(b.*FTYPI  (IM  ))-l  .) 
Eb=EXP(5.*FN)*<:>.*FN-l.  ) 

TC2».0l+.03b/<  l.-FTYPI(I«>) 

TC3=.007/(  »  ,-FTYPI  (  IM ) ) *(  E5-E'->F  )/(EXP(S  .  *FN  )-EXP(  b.  *FTYP1  (IM)  )) 
TOTCi  +TC2-TC3 

C  CALCULATE  CONSTANT  M  ,0R  ACTIVE  STATE  FNC  AND  MAX  STRAIN  RATE  C VEL) 

10  FM®.3/2/TC 

EMAX*.29//TC 

C  CALCULATE  PASSIVE  SPRING  FORCE 

FSP 31 . 00  1  6296* ( EXP (  / .  66  l  6*0)  - 1  .  ) 

C  CALCULATE  ACTIVE  STATE  FNC 

0=32. o3*( I .-EXP(-rM*TIME)) **2 . 

F0=(  .'J0b+0)/(  I  ,+Q) 

R=R  ATE/EMAX 

C  CALCULATE  LENGTH-TENSION  RELATIONSHIP 

FL=.32+.7I*EXP(-1  .  M2*D)*SIN<3.  /22* ( D+. 3444  )) 

C  CALCULATE  FORCE-VELOCITY  RELATIONSHIP 

FV=. 1 433/ ( . 10 /3+EXP ( - 1 . 409*5 IWIK 3.2 *R  +  I  .6 ))) 

C  CALCUL\TE  VISCOUS  DAMPING  FORCE 
FDAMP-.O 

I  F(  R.GT.O.  )FDAMP* .  1 .13*11 
C  CALCULATE  MAX  MUSCLE  FORCE 
FMAX  =  . 00673**( I  ,-FN) 

F*( FSP+FQ*FL*FV+FDAMP )*FMAX 
I F( F . LG . . 1 5 )G0T07b 
C  CALCULATE  ENDURANCE  TIME 

TEND*!  236 ,5/(  F*  I  Oj  .  -  I  b .  >**.  6 1  B  -  12. v 
I F ( *  END . LE . 6. 9  )TE!lD=£>  .9 
I F(  ri:lE.LE.rEND)Uoro7b 

C  If  TIME  EXCEEDS  ENDURANCE  TIME,  FORCE  CAPABILITY  iS  REDUCED 
T=(  TI  ME-TEND)  /60. 


B-l 


F=F*( 93 . 1 -23. 9*T+ 1 . 9*T*T) / 1 00. 

IF(F.LE.. I5)F».1d 

75  ELOSS  -  RATE*(FO*FL*FV)*FSCALE( IM)*F«AX 
0  GOLGI  TENDON  ORGAN  RESPONSE  (CLASH  KNIFE  REFLEX ) 

I F(  XTIME.EO.U.  uoToao 
DELr*TIME-TPREV( IM) 

IF(DFLT.EO.O.  )G0T098 
DELF*(F -FPREV ( I M )  )/DELT 

NRI  TE(  I  7, B9) TIME,  fPRGV(  IM) . DGLT. F, FPHGV  (IM)  ,DELF 
59  FORMAT!  6F  I O.ti) 

IF( 0ELF.UT.5o. )F=0. 

30  FPREV ( IM)*F 

TP.RGV(  I  M>**TIME 
GOTO 99 

9  9  F*FPREV  ( I M  ) 

99  FaF  *F SCALE ( Iw ) 

RETURN 

END 
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SOURCE  DECK  FOR 


con  o  c  c.  c.  o  o  c  o  o  o  c 


SUBROUTINE  H INPUT 


RE*/  6/3U/3~j 

CONTROLS  THE  INPUT  OF  CARDS  H.8.A  -  F.8.D  CONTAINING  THF  SETUP  AND 
CONTROL  OF  THE  HARNESS  BELT  SYSTEm. 


IMPLICIT  REAL*8(A-H,0-Z> 

COMMON/ CONTRL/  TlME.NSFG.NJNT .NPL.NBLT,  'AG.NVEH.NGRND, 

*  NS . NO , NSD.  NFLX .  NHRNSS,  Nrt  i  „DF ,  NJNTF , NPR  T ( 36  ) 
COMMOM/CNSNTS/  P I , RADIAN, G, TH I RD, EPS (24 ) , 

*  UN ITL, UN ITM .UNITT.GRAVTY! 3) 

COM MON/ HRN ESS/  BAR  ( I  b ,  I  00 ) ,  BB(  I  00) ,  BUDOT ( I  00 ) ,  PLOSS  ( 2  ,  I  00 ), 

*  XLONJ! bO ) ,HTI ME (2 ) , I BAR! b , I00),NL(2, IOj) , 

*  NPTSPB(bO) , NPTPLX ( bO ) , NTH R NS ( bO ) , NbLTPH ( o  J ) 

*  ,FSCALE( 50) .FTYPKjO) ,FMU(bO >, RECTI M(bU> ,QEL(bG) 
COMMON/TAtiLES/MXNTI  ,MXNTB,MXT3I  , MXTB2.NTI  (bO)  ,NTAb(  1230)  ,TAB(6000) 
C0MM0N/CNT5RF/  PL< I  7, 30), BELT ( 20,3) ,TPTS (6, 8 ) , BD(24 ,40 ) 

COMMON/TI TLES/  DATE ( 3) ,  CO.MENT (40 )  PSTTL! 20)  .uDYTTL (b ) , 

*  BLTTTL(b ,3 ) , PLTTL( 3,30), BAUTTL! b,6) , SEo! 30) , 

*  J0INT(30),CG3(30) ,JS(3J) 

REA L  D ATE . COMENT , VPSTTL , BDYTTL , BLTTTL, PLTTL , BAoTTL. SGw , JO  INT 
LOGICAL*!  CGS.JS 

THIS  COMMON/TEMTVS/  IS  SHARED  BY  CINPUT,  FINPUT,  HINpUT  AND  FDINIT 
COMMON/TEMPVS/  JTI TLEtb ,S I ) ,Nr( b) ,MS(3) , KTITLE ( 31 ) 

REAL  JTITLE.KTITLE 
IP  (NHRNSS. EO.O)  GO  TO  99 

INPUT  CARD  P.B.A 

(NOTEi  NHRNSS  NOW  SUPPLIED  ON  INPUT  CARD  D. I ) 

NbLTPH  -  NO.  OF  BELTS  PER  HARNESS 

READ< lb,  1 1  )  (NuLTPH(I) ,1-1 .NHRNSS) 

11  FORMAT! 1814) 

W R I TE ( 16,12)  NHRNSS, (NULTPH( I ), 1  =  1 , NHRNSS ) 

12  FORM AT ( '  I  HARNESS-BELT  SYSTEM  INPJT'.93X, 'CARDS  F.8'// 

*  '  No.  OF  HARNESSES  »',I4// 

*  '  NO.  OF  BELTS  PER  HARNESS  »',dOI2) 

Jl  »  I 
<1  =  1 

no  90  1=1, NHRNSS 
IF  (NBLTPH! I) .LE.O)  GO  TO  90 
J2  =  Jl  +  NBLTPH! I )  -I 

INPUT  CARD  F.8.B  -  NPTSPj  -  NO.  OF  POINTS  PER  BELT. 

READ  ( I  •_> ,  I  / )  (NPTSPB(J), J=JI ,J2) 

17  FORM AT( 1014) 

WRITE( I  6. 1 3)  I , (MPTSPB( J), J=JI ,J2  ) 

1.3  FORMAT! '0  FOR  HARNESS  NO.',  13, '  No.  Or  POINTS  PER  BELT  ='.2014) 

DO  30  J  =J I , Jp 

IF  ( NPTSPf)  ( J  ) .  EQ.O )  GO  TO  80 

INPUT  CARD  F.8.C  -  b  FUNCTION  NOS  AND  LENGTH  OF  EACH  BELT. 

READ! I b , 1 4 ) NF , XLONG! J ) , FSCALE ( J ) ,  rTYP I ( J ) , FMU! J ) , RECTI M( J ) , DEL! J ) 

I  4  j- ORM AT (  51 4 , F 1 2 . 0 .  bF3 . 0 ) 

C  CHEC<  IF  MUSCLE  CHARACTERISTICS  ARE  REASONABLE 
I F( RECT I M ( J ) . LT.  . I )RECTIM( J )= . I 
IF!  FTYPI (J).oO. I Oo.)FTYP» ( J)  =  99.9y 
IF! rMU! J) .LT.FTYPI (JJ)FMU! J ) =HT  Y  F I (J) 


IF!FMU! J) .LE.30.  .AND.  FMU<  J)  .GT.r-TYPI  ( J  ) )  FTYPI  (J)=i  mU!J  ) 

IF!  FTXPI (J).LE.O.O)  FTYPI (J)-.OOOi 

IF!  .“IF!  I  ) .GE.O) WRITE! I  o,  lb)  I .  J,ilr  ,X  LONG  (J  )  .UNITL 

15  FORMAT!  'O  HARNESS  NO.",  13,'  BELT  NO. '.13,'  FUNCTION  !IOS.'.:>I6. 

*  '  REFERENCE  SLACK  =  ',F9.3> 

IF!  NF(  I  >.LT.O>NRITE!  16.  19  )  I ,  J. FSCALE!  J ) , 

*1 J  N  l  TM ,  F  TY  P I  !  J  ) ,  F.M  J !  J  > .  R  ECT I M !  J ) ,  DEL  !  J ) 

19  FORMAT! 'O  MUSCLE  MO. '.1 3,'  SECTION  NO. '.1 3, 

*'  MUSCLE  FORCE  SCALING  FACTOR«',r3.2, 1  X, A4/ 

*'  3  SLOW  TWITCH  MU*',Fd.2,'  %  MU  RECRUI TED»',Fd.2, 

*'  RECRUITMENT  TIME®'  ,FU.4,'  DELAlf  TIME=',F0.4) 

IF  ! XLONG! J ) . EO.U ,U )  XLONoIJ)  -  EPST24) 
rTYPI ! J )-FTYPI ! J )/! 00. 

FM’J !  J)=FMU!  J) /  I  00. 

WRITE!  16,  16) 

16  FORMAT  !'U  K  KS  KE  NT  NPD  NDr  FUNCTION  NOS. 

*  66 X,' CARDS  F.B.D'/) 

SET  UP  POINTERS  IN  NTAd  AND  INITIAL  VALUES  OF  TAB  FOR  BELT  J 
AS  WAS  DONE  FOR  OTHER  CONTACTS  IN  SUBROUTINE  FINi'UT. 

NTHRNS! J)  «  MXNTb+l 
CALL  FDINIT 

K2  -  rJ  ♦  NPTSPB! J )  -  I 
DO  10  K-KI.K2 

INPUT  CARD  F.S.D 

READ! I  5 ,2  I  )  KS ,  .«E ,  NPD ,  NDR , MF ,  (UAR!L,K  ),L-1 ,3) 

21  FORMAT  ! 9 1 4 » 3F 1 2 . U ) 

READ! 15,22)  ! BAR! L ,K ) , L-V , »2 ) 

22  FORMAT  (6FI2.U) 

IBAR! I ,K)  «  KS 

I B A R ! 2 ,  K )  «  KE 
I  BAR! 4, K)  ®  NPD 
IRAR! 5, K)  =  NDR 
IBAR!3,K)  -  MXNTB+I 
CALL  FDINIT 
SORER  »  I  .0 

IF  IkE.NE.O)  SORER  »  DSQRT ! XDY! BAR! I,  K ) « BD! / . KE ) . BAR !  I  , K >  > ) 
DO  26  L-1,3 

IF  (KE.NE.O)  SAR!L+6,K)  -  BO! L+3, KE) 

26  BAR ( L+3.K )  »  UAR(L,K)/SORER 

WRITE! 16,31 )  K, !  IiiAE!L,<  >  ,L-i .  5>  ,NF 
31  FORMAT  (I  116) 

70  CONTINUE 

WRITE!  I  6,  VI  )  UNITL, IJNITL.UNITL, UNITL 

71  FORMAT  !'0', I 2X, 'BASE  REFERENCE  !',  A4,')'. 

*  !X, 'ADJUSTED  REFERENCE  A4,')', 

*  IIX, 'OFFSET  !',  A4, ') ' , 

*  I  I  X, 'PREFERRED  DIRECTION  !',A4,')'/ 

*  5X,'K',  4!8X,'X',0X.'Y',8X, 'Z',3X)  /) 

WRITE! 16, 12)  (K, (uAR(L.K) ,L-l , I2),K-KI ,K2) 

72  FORMAT  !  1 6,  3X. 3Fv> . 3 , 3X.  3F9. .1 , 3X , 3F9 . 3,  3X .  3F9. 3 ) 

K  I  =  ;» 2  ♦  I 

80  CONTINUE 
Jl  *  J2  +  I 
90  CONTINUE 

DO  92  K-1,100 

BBDOTIK)  -0.0 


APPENDIX  D 


VALIDATION  OF 
VOL. 


I 

i 


REVISIONS  TO 

THE  CRASH  VICTIM  SIMULATOR 
3  -  USER'S  MANUAL 


APPENDIX  D 


REVISIONS  TO 

VALIDATION  OF  THE  CRASH  VICTIM  SIMULATOR 
VOL.  3  ~  USER'S  MANUAL 

D.  SUBROUTINE  SINPUT  (p.  46) 

NHRNSS  Number  of  harness-belt  systems  or  active  muscles  to  be 

supplied  on  cards  F.8.B-F.8.D.  May  be  zero  or  blank. 
Maximum  value  -  50.  Note:  In  version  12  (for  WPAFB) 
this  variable  was  supplied  on  card  F.8.A. 

F.8  Subroutine  HINPUT  -  card  input  for  harness-belt  systems  or  active 

musculature  (p.  69) 

Note:  NHRNSS  which  was  supplied  on  card  F.8. A  for  version  12  is  now  supplied  on 

card  D. 1 .  If  NHRNSS  #0,  cards  F.8  must  be  supplied.  Previously  for  version  12, 

a  blank  card  F.8. A  was  required  if  no  harness  belt  systems  were  desired. 

CARD  F.8. A  FORMAT  (5014)  use  three  cards  if  NHRNSS  >  18 

NBLTPH(I),  Number  of  individual  belts  for  each  harness  No.  1.  May 

1-1 , NHRNSS  be  zero  or  blank.  Maximum  value  of  sum  of  all  NBLTPH  is 

50. 

Card  F.8. A  is  followed  by  NHRNSS  sets  of  cards  F.8.B  -  F.8.D. 

Card  F.8.B  Format  (1814)  use  two  cards  if  NBLTPH(I)  *  18. 

NPTSPB(J),  The  number  of  reference  points  including  anchor  points 

J»1 ,NBLTPH( I)  for  belt  No.  J  of  harness  No.  1  may  be  zero  or  blank. 

The  maximum  value  of  the  sum  of  all  NPTSPB  for  all 
harness-belt  systems  is  100.  The  maximum  value  of  the 
sum  of  all  NPTSPB  for  any  one  harness  belt  system  is  50. 
The  maximum  value  of  any  individual  NPTSPB  is  25. 

Each  card  F.8.B  is  followed  by  NBLTPH(I)  sets  of  cards  F.8.C  -  F.8.D 

Card  F.8.C  Format  (514,  F12.0,  5F8.0) 

NF(L),L«1  The  function  numbers  from  cards  E.1  to  define  the 

stress-strain  of  belt  No.  J.  The  definition  of  these 
functions  are  identical  to  those  of  NF(1)  to  NF(5)  on 
cards  F.2.B,  except  that  the  use  of  rate  dependent 
functions  is  permitted.  NF(1)  positive  indicates  a 
harness  while  NF(1)  negative  indicates  active 
musculature.  The  function  numbers  are  not  used  for  the 
active  musculature. 

The  initial  slack  (IN)  of  belt  No.  J.  A  negative  value 
can  be  specified  to  indicate  a  pre-tightened  belt.  The 
program  will  add  this  to  the  initial  geometric  length  to 
obtain  the  initial  belt  length  and  distribute  the  slack 


XI  0NG(J) 


FSCALE(J) 


FTYP1 ( J) 


FMU(J) 


RECTIM(J) 


DEL(J) 


proportionately  between  the  points.  For  the  active 
musculature  the  slack  is  not  needed  (i.e.,  negative 
strain  is  allowed)  and  can  be  left  blank. 

Scaling  used  to  multiply  the  normalized  forces  of  each 
muscle,  equivalent  to  the  maximum  force  the  muscle  can 
produce. 

Percent  of  Type  1  motor  units  or  muscle  fibers  in  each 
muscle.  The  value  of  FTYP1  must  be  equal  to  FMU  for  FMU 
<  30 S  and  less  than  or  equal  to  FMU  if  FMU  >  30%. 

Percent  of  motor  units  recruited.  FMU  can  vary  from  0%, 
or  no  muscle  activity,  to  1 00%  or  maximum  voluntary 
contraction. 

Time  for  full  motor  unit  recruitment  in  seconds.  RECTIM 
must  be  greater  than  .1  sec  found  in  ballistic 
contractions. 

Delay  time  before  motor  unit  recruitment  in  seconds. 


Table  7.  ATB  Specifications  for  Shoulder  Musculature  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (in)  Scaling 

Factor 


Muscle  Group  1 

SEC 

X 

Y 

Z 

SEG 

X 

Y 

Z 

(lbs) 

1) 

Deltoid  a) 

3 

-.79 

6.49 

-3.81 

12 

.39 

.39 

.07 

121.11 

b) 

3 

.79 

6.49 

-3.81 

12 

.39 

.39 

.07 

121.11 

2) 

Supraspinatus 

3 

.0 

3-74 

-1 .85 

12 

.0 

1.57 

-5.83 

72.6 

3) 

Pectoralis  a) 

3 

2.0 

.0 

-.28 

12 

•  39 

•  39 

-1.9 

74.8 

Major  b) 

3 

2.0 

3.34 

-2.64 

12 

•  39 

•  39 

-1.9 

74.8 

4) 

Latisslmus  a) 

2 

-2.0 

.0 

3.85 

12 

.0 

.39 

-2.68 

59.18 

Dorsi  b) 

3 

-2.0 

.0 

.71 

12 

.0 

.39 

-2.68 

59.18 

5) 

Teres  Major 

3 

.0 

4.92 

2.09 

12 

.0 

-.39 

-3.08 

109.34 

6) 

Teres  Minor 

3 

.0 

5.71 

.51 

12 

.0 

.39 

-5.44 

34.54 

7) 

Subscapularis 

3 

.39 

4.92 

.51 

12 

.59 

.59 

-5.44 

217.8 

8) 

Coraco  Brachialis 

3 

1.57 

7.28 

-3.23 

12 

.0 

-.39 

-.72 

33.44 

9) 

Infraspinatus 

3 

-.39 

4.52 

-.28 

12 

.0 

.59 

-5.44 

131 .56 

Table  8.  ATB  Specifications  for  Hip  and  Knee  Musculature  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (in)  Scaling 

Factor 


Muscle  Group 

SEG 

X 

Y 

Z 

SEG 

X 

Y 

Z 

(lbs) 

1 )  Gluteus  Medius 

1 

-1.46 

.0 

-2.23 

6 

-1.97 

5.88 

-8.21 

465.96 

2)  Gluteus  Minimus 

1 

.0 

.0 

-1.44 

6 

-1.18 

5.88 

-8.21 

21 1 .2 

3)  Tensor  Fasciae  Latae 

i 

-1 .57 

.0 

-2.23 

6 

-1.97 

5.49 

-5.26 

54.56 

4)  Obturatus  Internus 

1 

.0 

1.95 

3.68 

6 

-1.18 

3.49 

-8.88 

86.02 

5)  Adductor  Longus 

1 

.0 

.76 

3.28 

6 

-1.18 

4.7 

-1.32 

110.66 

6)  Adductor  Brevis 

1 

.0 

1.55 

3.68 

6 

-1.18 

4.7 

-3-09 

99.88 

7)  Adductor  Magnus  a) 

1 

.0 

3-13 

3.68 

6 

-1.18 

4.7 

-6.44 

226.38 

b) 

1 

.0 

1.55 

3.68 

6 

-1.18 

4.7 

-8.8 

226.38 

8)  Pectineus 

1 

.787 

1.55 

.92 

6 

-1.18 

1 .0 

-5.06 

54.34 

9)  Quadratus  Femoris 

1 

-.79 

1.95 

3.68 

6 

-2.36 

-1.97 

-7.62 

64.02 

10)  Obturatus  Externus 

1 

.0 

1.95 

3.68 

6 

-1.18 

-1.97 

-8.8 

108.9 

11)  Gluteus  Maximus  a) 

1 

.0 

-3-15 

-.65 

6 

-1.18 

1.38 

-4.47 

323.62 

b) 

1 

.0 

.0 

.0 

6 

-1.18 

1.38 

-4.47 

323-62 

12)  Semimembranosus 

1 

-1.57 

.0 

2.49 

.0 

-1.18 

-7.14 

285.34 

13)  Semitendinosus 

1 

-1 .97 

.0 

2.49 

.39 

-.79 

-4.39 

95.26 

14)  Biceps  Femoris  a) 

-1.77 

1.18 

2.49 

.0 

1  .57 

-7.14 

129.8 

b) 

-1.18 

1.18 

-.53 

.0 

1  .57 

-7.14 

129.8 

15)  Quadriceps  a) 

.0 

.0 

-.23 

1.18 

.0 

-7.93 

308.0 

b) 

.0 

1.18 

.65 

1.18 

.0 

-7.93 

924.0 

16)  Iliopsoas  a) 

1 

.0 

-.39 

3.K1 

-1.18 

.79 

-6.44 

165.66 

b) 

2 

.0 

.0 

-.17 

-1.18 

.79 

-6.44 

165.66 

17)  Gastrocnemius 

6 

.0 

.0 

8.52 

8 

.0 

.0 

-1  .5 

344.52 

18)  Popliteus 

6 

-1.18 

-1.57 

9.7 

7 

-.79 

.0 

-5.57 

43.78 

19)  Gracilis 

1 

-.39 

.76 

3.67 

7 

.39 

-.79 

-5.17 

35.86 

20)  Sartorius 

1 

-1.57 

.0 

-1 .05 

7 

.39 

-.79 

-5.17 

34.1 

i 


i 

3 
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Table  9  ATB  Specifications  for  Neck  Musculature  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (in)  Scaling 


Factor 


Muscle  Group 

SEC 

X 

Y 

Z 

SEG 

X 

Y 

z 

(lbs) 

1) 

Longus  Capitis 

4 

-.75 

1.0 

-.54 

5 

1.54 

.39 

3.18 

16.5 

2) 

Rectus  Capitis 
Anterior 

4 

-.75 

1.0 

-2.52 

5 

1.26 

.39 

3.18 

5.5 

3) 

Rectus  Capitis 
Lateralis 

4  | 

-.75  | 

1  .0 

-2.52  | 

5  I 

.87  | 

1 .18| 

3.18 

5.5 

4) 

Rectus  Capitis 
Posterior  Major 

4  j 

-.75  | 

.0 

-1.86 

5 

-.55  | 

.94| 

2.83  | 

11.0 

5) 

Rectus  Capitis 
Posterior  Minor 

4 

-.75  | 

.0  I 

-2.52  j 

5 

-.55  | 

.  3l| 

2.98 

;  8.47 

6) 

Obliquus  capitis 
Superior 

4 

I  ”*75 

1.0  I 

-2.52 

'  5 

^  5 

-.55 

1.38 

2.67 

j  22.0 

7) 

Splenius  Capitis 

4 

1 

.0 

2.16 

-.39 

1 .57 

2.67 

26.84 

8) 

Longisslmus  Capitis 

4 

-.75 

1.0 

2.16 

5 

.08 

1 .34 

2.67 

1 1.0 

9) 

Spinalis  Capitis 

4 

-.75 

1.C 

2.16 

5 

-.87 

.24 

2.83 

1 1.0 

10) 

Semispinalis  Capitis 

4 

-.75 

1.0 

2.16 

5 

1  -1.42 

.79 

2.67 

52.36 

11  ) 

Trapezius 

3 

.0 

7.28 

-2.35 

5 

-1.42 

.39 

2.59 

233.2 

12) 

Sternocleido- 

Mastoideus 

3 

3.0 

.0 

-2.0 

5 

■  ".39 

|  1.57 

2.67 

35.2 

13) 

Levator  Scapulae 

3 

-2.0 

4.28 

-1.85 

2 

-.31 

.45 

-1.85 

390.5 

14) 

Longus  Colli 

3 

.75 

.0 

-2.0 

-  4 

.75 

.0 

-.18 

16.5 

15) 

Scalenus 

3 

3.0 

2.0 

-1.85 

4 

-.75 

1  .0 

.2 

38.5 

16) 

Splenius  Servicis 

3 

-.75 

.0 

-1.26 

4 

-.75 

1  .0 

-1.86 

15.4 

17) 

Longisslmus  Servicis 

3 

-.75 

1 .0 

-1.26 

4 

-.75 

1.0 

-.54 

13.2 

18) 

Spinalis  Cervicis 

3 

-.75 

.0 

-2.6 

4 

-.75 

.0 

-1.86 

27.5 

19) 

SemiSpinalis 

Cervicis 

3 

-.75 

,.o 

-2.6 

4 

-.75 

l  .0 

I 

-1.2 

44.0 
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Table  10  ATB  Specifications  for  Trunk  Musculature  (Right  Side  Only) 


ORIGIN  (in)  INSERTION  (in) 


Muscle  Group  1 

SEd 

X 

Y 

Z 

SEG 

X 

Y 

z 

1) 

Illocos tails 

Lumborum 

1 

-1.0 

1.5 

3.85 

2 

-1.0 

3.0  | 

.31 

2) 

Iliocostalis  Dorsi 

2 

-1.0 

.31 

3 

-1.0 

3.0 

.0 

3) 

Longisslmus  Dorsi 

2 

-1.0 

-.19 

3 

-.75 

1 .0 

.71 

4) 

Spinalis  Dorsi 

2 

-1.0 

-4.23 

3 

-.75 

1 .0 

.71 

5) 

Semlsplnalls  Dorsi 

3 

-.75 

1.89 

4 

-.75 

1.0 

1.5 

6) 

Multifldus  a) 

2 

-1.0 

-2.89 

3 

-.75 

0 

4.23 

b) 

3 

-.75 

-2.04 

4 

-.75 

0 

2.82 

7) 

Intersplnales  a) 

2 

1.0 

.0 

-2.89 

3 

-.75 

0 

4.23 

b) 

3 

-.75 

.0 

-2.04 

4 

-.75 

0 

2.82 

8) 

Iniertrans-  a) 

Versarii 

2 

-1.0 

1.5 

-2.89 

3 

-.75 

1.0 

4.23 

b) 

3 

-.75 

1.0 

-2.04 

4 

-.75 

1.0 

2.82 

9) 

External  Oblique 

1 

3.0 

3.52 

-3.08 

2 

.0 

3.52 

.31 

10) 

Internal  Oblique 

1 

.0 

3.52 

-3.08 

2 

3.0 

3.52 

.31 

11  ) 

Rectus  Abdominus 

1 

1.0 

.5 

1.92 

3 

3.0 

.5 

-4.23 

12) 

Quadratus  Lurborum 

1 

.0 

3.52 

-3.08 

2 

-.75 

1.0 

-.17 

Scaling 

Factor 

(lbs) 


22.0 

11.0 

22.0 

22.0 

22.0 

12.75 

13-75 

5.5 

5.5 

2.75 

2.75 

150.7 

12*1.95 

58.52 

61.6 
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